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The existence of dark matter was suggested, using simple gravitational arguments, sev- 
enty years ago. Although we are now convinced that most of the mass in the Universe is 
indeed some non-luminous matter, we still do not know its composition. The problem of 
the dark matter in the Universe is reviewed here. Particle candidates for dark matter are 
discussed with particular emphasis on Weakly Interacting Massive Particles (WIMPs). 
Experiments searching for these relic particles, carried out by many groups around the 
world, are also reviewed, paying special attention to their direct detection by observing 
the elastic scattering on target nuclei through nuclear recoils. Finally, we concentrate on 
the theoretical models predicting WIMPs, and in particular on supcrsymmctric exten- 
sions of the standard model, where the leading candidate for WIMP, the neutralino, is 
present. There, we compute the cross section for the direct detection of neutralinos, and 
compare it with the sensitivity of detectors. We mainly discuss supergravity, superstring 
and M-theory scenarios. 
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1. Introduction 

One of the most evasive and fascinating enigmas in physics is the problem of the 
dark matter in the Universe. Substantial evidences exist suggesting that at least 
90% of the mass of the Universe is due to some non-luminous matter, the so-called 
'dark matter'. However, although the existence of dark matter was suggested 70 
years ago, we still do not know its composition. 

In 1933, the astronomer Fritz Zwicky ^ provided evidence that the mass of 
the luminous matter (stars) in the Coma cluster, which consists of about 1000 
galaxies, was much smaller than its total mass implied by the motion of cluster 
member galaxies. However, only in the 1970s the existence of dark matter began 
to be considered seriously. Its presence in spiral galaxies was the most plausible 
explanation for the anomalous rotation curves of these galaxies, as we will discuss 
in detail in Section 2. In summary, the measured rotation velocity of isolated stars 
or gas clouds in the outer parts of galaxies was not as one should expect from the 
gravitational attraction due to the luminous matter. This led astronomers to assume 
that there was dark matter in and around galaxies. We will see that nowadays 
there is overwhelming observational evidence for the presence of dark matter. It 
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not only clusters with stellar matter forming the galactic halos, but also exists 
as a background density over the entire Universe. Thus the problem is no longer 
to explain the rotation curves but to decipher the nature of the dark matter. In 
any case, it is fair to say that some authors still suggest that dark matter is not 
necessary to explain the rotation curves. They try to avoid the introduction of this 
hypothesis modifying the Newton's law at galactic scales. Although these attempts 
are not very convincing, we will briefly discuss them 

As we will explain in detail in Section 3, the search of the dark matter provides a 
potentially important interaction between particle physics and cosmology, since only 
elementary particles are reliable candidates for the dark matter in the Universe. In 
particular we will see that baryonic objects, such as e.g. gas, brown dwarfs, etc., can 
be components of the dark matter, but more candidates are needed. The reason 
being that they cannot be present in the right amount to explain the observed 
matter density of the Universe, 0.1 < flh 2 < 0.3. Fortunately, particle physics, and 
mainly extensions of the standard model offer candidates for dark matter. Indeed, 
detecting non-baryonic dark matter in the Universe might be a signal for new 
physics beyond the Standard Model. It is fair to say that many of these candidates 
are quite exotic and most of them are ephemeral. However, we will also see that 
very interesting and plausible candidates for dark matter are Weakly Interacting 
Massive Particles (WIMPs), since long-lived or stable WIMPs can be left over from 
the Big Bang in sufficient number to account for a significant fraction of relic matter 
density. As suggested in 1985 by Goodman and Witten^, and also by Wasserman 
121 this raises the hope of detecting relic WIMPs directly. 

WIMPs cluster gravitationally with ordinary stars in galactic halos, and there- 
fore they must be present in our own galaxy, the Milky Way. As a consequence 
there will be a flux of these dark matter particles on the Earth, and experiments 
can be carried out to detect this flux by observing its elastic scattering on target 
nuclei through nuclear recoils. As we will mention, inelastic scattering is also possi- 
ble but the event rates seem to be suppressed. Thus in Section 4 we will review the 
current and projected experiments for detecting WIMPs directly in underground 
labora torie s. In fact, a dark matter experiment of this type, the DAMA collabo- 
ration EE1 has reported data recently favouring the existence of a WIMP signal. 
It was claimed that the preferred range of the WIMP-nucleon cross section is rs 
10~ 6 - 10~ 5 pb for a WIMP mass of about 100 GeV. Although this could be a great 
discovery, it is fair to say that this experiment is controversial, mainly because other 
collaborations reported negative search results. In particular, CDMS ^, and more 
recently EDELWEISS^ and ZEPLIN 1^, claim to have excluded important regions 
of the DAMA parameter space. In any case, the whole region allowed by DAMA 
will be covered by these and other experiments in the near future. Besides, there 
are projected experiments which will be able to measure a cross section as small as 
10~ 9 pb. In addition, there are also promising methods for the indirect detection of 
WIMPs by looking for evidence of their annihilation, as we will briefly discuss. Thus 
it seems plausible that the dark matter will be found in the near future. Given this 
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situation, it is crucial to analyze in detail the present theoretical models predicting 
WIMPs. 

In Section 5 we will review the situation concerning this point. In particular, we 
will concentrate on the leading candidate for WIMP, the li ghtes t supersymmetric 
particle (LSP). As is well known, supersymmetry (SUSY) lallUI is a new type of 
symmetry, since it relates bosons and fermions, and this ensures the stability of 
the hierarchy between the weak and the Planck scales. This is the most relevant 
theoretical argument in its favour. In SUSY models the so-called R-parity is of- 
ten imposed in order to avoid fast proton decay or lepton number violation. This 
yields important phenomenological implications. SUSY particles are produced or 
destroyed only in pairs and therefore the LSP is absolutely stable, implying that 
it might constitute a possible candidate for dark matter. So SUSY, whose original 
motivation has nothing to do with the dark matter problem, migh t so lve it. 

The first discussion of SUSY dark matter was by Goldberg El m 1933. He 
considered the SUSY fermionic partner of the photon, the photino, as the LSP 
and pointed out the strong constraints on its mass from its relic abundance. Soon 
followed two works considering also this possibility l -^ l -^ l but it was in Ref. fPTl 
were the analysis of the most general case was carried out. In particular, one has 
to take into account that in SUSY the photino mixes with the fermionic partners 
of the Z° and the two neutral Higgs bosons. Therefore one has four particles called 
neutralinos, x° (i = 1,2,3,4). Obviously, the lightest neutralino, Xit will be the 
dark matter candidate. The fact that the LSP is chosen to be an electrically neutral 
particle (also with no strong interactions) is welcome since otherwise it would bind 
to nuclei and would be excluded as a candidate for dark matter from unsuccessful 
searches for exotic heavy isotopes. It is remarkable that in large regions of the 
parameter space of the simplest SUSY extension of the standard model, the so- 
called Minimal Supersymmetric Standard Model (MSSM), the LSP is indeed the 
lightest neutralino. 

It is worth noticing here that several accelerator experiments are in preparation 
in order to detect SUSY particles. For example, LHC at CERN will probably start 
operations in 2007 producing energies about 1 TeV. Concerning the LSP, one will 
be able to detect events with SUSY particle decays which produce lots of missing 
energy. However, let us remark that, even if neutralinos are discovered in this way, 
only their direct detection in underground experiments due to their presence in our 
galactic halo will confirm that they are the sought-after dark matter of the Universe. 
As a matter of fact, given the present and projected experiments in underground 
laboratories in order to detect WIMPs, there may be a competition between these 
experiments and those in accelerators in the hunt for the first SUSY particle. 

Thus we will review how big the cross section for the direct detection of neu- 
tralinos can be, in the generic context of the MSSM. In fact, the cross section for 
the elastic scattering of x® o n nu cleons has been examined exhaustively in the lit- 
erature since many years ago^l Obviously, this computation is crucial to analyze 
the compatibility of the neutralino as a dark matter candidate, with the sensitivity 
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of detectors. In the light of the recent experimental results mentioned above, many 
theoretical works studying possible values of the cross section in different scenarios 
have appeared. Here we will concentrate on the mo s t re cent results in supergravity 
(SUGRA), superstring and M-Theory scenarios UfiHSSI 

Let us recall that in the framework of SUGRA * the scalar and gaug- 
ino masses and the bilinear and trilinear couplings are generated once 
SUSY is softly broken through gravitational interactions, and in addi- 
tion radiative electroweak symmetry breaking is imposed. First, we will 
analyse this framework in the usual context of a Grand Unified The- 

ory (GUT) with scale Mgnr ~ 2 x 10 16 GeV and universal soft terms 

18 19 21 22 23 24 26 28 29 31 36 38 39 40 43 44 45 46 47 48 49 50 51 55 57 58 59 60 63 64 65 

* 

the so-called minimal supergravity (mSUGRA) scenario. Second, we will discuss 
how the results are modified when the GUT condition is relaxed. In particular, 
we will consider the case of an intermediate scale | 33 | 35 | 41 | 45 | 5U | 63 ] a possibility 
inspired by experimental observations and also by some superstring constructions. 
Then, a more general situation in the context of SUGRA, non-universal scalar 

\im\\m±mmW\Am)\immm® and gaugin0 mum^ssm^m masseS) 

will be studied. On the other hand, a phenomenological SUSY model, with pa- 
rameters defined directly at the ele ctroweak scale, the so-called effMSSM, has also 
been studied in the literature We will discuss this approach briefly. 

Let us remark that in these analyses to reproduce the correct phenomenology is 
crucial. In this sense it is important to check that the parameter space of the dif- 
ferent scenarios fulfills the most recent experimental and as trop hysical constraints. 
Concerning the former, the lower bound on the Higgs massif, the b — > S7 branch- 
ing ratio and the muonEHI g — 2 will be taken into account. The astrophysical 
bounds on the dark ma tter density mentioned above will also be imposed on the 
theoretical computation of the relic neutralino density, assuming thermal pro- 
duction. In addition, the constraints that the absence of dangerous charge and 
colour breaking minima imposes on the parameter space will also be taken into 
account. 

Finally, since generically the low-energy limit of superstring theory i s 4_ 
dimcnsional SUGRA, the neutralino is also a candidate for dark matter in super- 
string constructions. Taking into account that the soft terms can in principle be 
computed in these constructions, we will review the associated Xi-nucleon cross sec- 
tion. In particular w e will discuss this first in the context of the (weakly-coupled) 
heterotic superstri ng ^ Then, we w ill consider D-brane configurations from the 
type I superstring ElEOEZEBBIBBEOI which are explicit scenarios where interme- 
diate scales, and also non-universality of soft terms, may occur. 

On the other hand . 4- dimensional SUGRA is also the low-energy limit of 11- 
dimcnsional M-Theoryl^J. The latter has been proposed as the fundamental theory 
which contains the five 10-dimensional superstring theories, and, in particular, com- 
pactified on an orbifold is indeed the strong coupling limit of the E$ x E% heterotic 
superstring theory. Therefore, the same review for the cross section as above will be 
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carried o ut for this strongly-coupled heterotic superstring obtained from M-Theory 
I ^U | 42 | 52 | -j_ ne g0 _ ca n ec j heterotic M-Theory. The conclusions are left for Section 6. 

2. The dark matter problem 

To compute the rotation velocity of stars or hydrogen clouds located far away 
from galactic centres is easy. One only needs to extrapolate the Newton's law, 
which works outstandingly well for nearby astronomical phenomena, to galactic 
distances. For example, for an average distance r of a planet from the center of the 
Sun, Newton's law implies that v 2 (r)/r = GM(r)/r 2 , where v(r) is the average 
orbital velocity of the planet, G is the Newton's constant and M(r) is the total 
mass inside the orbit. Therefore one obtains 



Clearly, v(r) decreases with increasing radius since M(r) is constant and given by 
the solar mass. 

In the case of a galaxy, if its mass distribution can be approximated as spherical 
or ellipsoidal, Eq. can also be used as an estimate. Thus if the mass of the galaxy 
is concentrated in its visible part, one would expect v(r) = y/ GM V i S / r oc l/^fr 
for distances far beyond the visible radius. Instead, astronomers, by means of the 
Doppler effect, observe that the velocity rises towards a constant value v c w 100 
to 200 km s _1 . An example of this can be seen in Fig. ^ (from Ref. I72J) . where 
the rotation curve of M33, one of the about 45 galaxies which form our small 
cluster, the Local Group, is shown. For comparison, the expected velocity from 
luminous disk is also show n. Thisp henomenon has already been observed for about 
a thousand spiral and in particular also for our galaxy, the Milky 

Way. Although this observation is more problematic in galaxies other than spirals, 
such as ellipticals, dwarf irregulars, dwarf spheroidals, lenticulars, etc., they also 
produce similar results ESEZIZSI 

The most common explanation for these flat rotation curves is to assume that 
disk galaxies are immersed in extended dark matter halos. Thus for large distances 
M{r)/r is generically constant because the mass interior to r increases linearly 
with r, M to t(r) = G~ 1 v 2 r. In fact, a self-gravitating ball of ideal gas at an uniform 
temperature of kT = ^m DM v 2 , where m DM is the mass of one dark matter particle, 
would have this mass profile In addition, dark matter is necessary to explain 
structure formation, as deduced from the cosmic microwave background measure- 
ments. Clumps of neutral particles arose first through gravitational attraction and 
later, when the neutral atoms were formed, these were gravitationally attracted by 
the dark matter to form the galaxies. 

The above analysis of rotation curves implies that 90% of the mass in galaxies 
is dark. Whereas current observations of luminous matter in galaxies determine 
n ium < 0.01, analyses of rotation curves imply in fact as 0.1. Let us recall that 
Q = p/pc is the present-day mass density averaged over the Universe, p, in units 




(1) 
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Fig. 1. Observed rotation curve of the nearby dwarf spiral galaxy M33, superimposed on its 
optical image. 

of the critical density, p c — 10 _5 ft. 2 GeV cm -3 , with h ~ 0.7. In fact, the previous 
value, f2 « 0.1, is really a lower bound, since almost all rotation curves remain flat 
out to the largest values of r where one can still find objects orbiting galaxies. We 
do not really know how much further the dark matter halos of these galaxies extend 
(see e.g. Fig.f^i. Thus we can conclude that galactic rotation curves imply > 0.1. 

Although the analysis of dark matter in cluster of galaxies become s more in- 
volved than in galaxies, different techniques have been used to compute the 
value of O. For example, the conventional method of studying the motion of cluster 
member galaxies seems to point at values of f2 larger than those obtained in galactic 
scales. The high temperature of the gas detected in clusters through its X-ray emis- 
sion also points at large amounts of dark matter. Finally, the more reliable method 
of studying the gravitational lensing confirms the previous conclusions. Here a clus- 
ter acts as a lens which distorts the light emitted by quasars and field galaxies in 
its background, due to the gravitational bending of light. All these analyses favour 
a value f2 w 0.2 — 0.3. Moreover, measurements of velocity flows at scales larger 
than that of clusters favour^ also very large amounts of dark matter: > 0.3. 
Thus the following astrophysical bounds are commonly used in the literature: 

0.1 < n DM h 2 < 0.3 . (2) 

As a matter of fact, the lower part of this rang e seems to be preferred by comparing 
models with a set of cosmological observations I ^ I ^ J. In particular, the (cold) dark 
matter range 

0.094 < fl DM h 2 < 0.129 , (3) 

can be deduced from the recent data obtained by the WMAP satellite 

Clearly, any sensible candidate for dark matter must be able to reproduce the 
observations summarized by the above equations, i.e. it must be present in the right 
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amount to explain the observed matter density of the Universe. This is the issue 
that we will review in the next section. 

Before finishing this section, it is fair to say that a small number of authors 
suggest that dark matter is not really necessary to explain rotation curves \U\ 
Basically their approach consists of modifying the Newton's law at galactic scales. 

OO OO OA 

This has been done first in the pure classical theory ' ° | OJ | ° ' For example, in 
Ref. I83lthe authors assume that the behaviour of the potential at galactic distances 
is V ~ lnr. In the context of the relativistic theory, several authors have s tudied 
the modification induced to the scalar potential by higher-derivative gravity ESESI 
or by a cosmological constant E3. These modifications have also been studied in 
other contexts, such as fiv e di mensional gravity inspired in strings ^Hl (D-branes), 
and geodetic brane gravityEHI. However, these attempts are not only rather ad hoc 
in general (the authors impose specific values for the free parameters of the theory 
in order to reproduce some of the rotation curves that have been observed) but 
also insufficient to account for the necessity of dark matter in scales larger than 
galactic ones^l (values of the parameters necessary to reproduce galactic rotation 
curves cannot reproduce the observations at larger scales). Recently, the authors 
of Ref. |^ using the Sloan Digital Sky Survey, have shown the first observational 
evidence that the halo density decline as 1/r 3 , as predicted by cold dark matter 
cosmological models. Alternative theories of gravity are in contradiction with this 
result. 



3. Dark matter candidates 

Here we will review possible candidates for dark matter. Since ordinary matter is 
baryonic, the most straightforward possibility is to assume also this composition for 
the evasive dark matter. The contribution from gas is not enough, so astrophysical 
bodies collectivelyknown as MAssive Compact Halo Objects (MACHOs) are the 
main candidates"^. This is the case of brown and white dwarfs, Jupiter- like ob- 
jects, neutron stars, stellar black hole remnants. However, the scenario of Big-Bang 
nucleosynthesis, which explains the origin of the elements after the Big Bang, taking 
into account measured abundances of helium, deuterium and lithium sets a limit 
to the number of baryons that can exist in the Universe, namely f2 baryon /i 2 < 0.05. 
This density is clearly small to account for the whole dark matter in the Universe 
(see bounds l|2*|l). The conclusion is that baryonic objects are likely components of 
the dark matter but more candidates are needed. This result is also confirmed by 
observations of MACHOs in our galactic halo through their gravitational lensing 
effect on the light emitted by stars. Their contribution to the dark matter density 
is small. Thus non-baryonic matter is required in the Universe. 

Particle physics provides this type of candidate for dark matter. In principle, the 
three most promising are 'axions', 'neutrinos' and 'neutralinos' with masses of the 
order of 10~ 5 eV, 30 eV and 100 GeV, respectively. As a matter of fact, neutrinos 
are the only candidates which are known to exist. However, although the other 
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particles are not present in the standard model, they are crucial to solve important 
theoretical problems of this model. Thus they are well motivated by extensions of 
the standard model. 

Before analyzing these three (standard) candidates, it is fair to say that many 
others have also been proposed. Although some of them are quite exotic, we will 
also discuss them in some detail below. 



3.1. Standard candidates 
3.1.1. Neutrinos 

In recent years, observation of solar and atmospheric neutrinos have indicated that 
one flavour might change to another EH. As is well known, these neutrino oscil- 
lations can only happen if neutrinos have mass. The best evidence for neutrino 
masses comes from the SuperKamiokande experiment in Japan concerning atmo- 
spheric neutrino oscillations. The results of this experiment indicate a mass differ- 
ence of order 0.05 eV between the muon neutrino and the tau neutrino. If there 
is a hierarchy among the neutrino masses (as it is actually the case not only for 
quarks but also for electron-like leptons), then such a small mass difference implies 
that the neutrino masses themselves lie well below 1 eV. This is not cosmologically 
significant as we will show below, since a light (m v < 100 eV) neutrino has a cos- 
mological density!^ f2„ ss m„/30 eV. On the other hand, there could be near mass 
degeneracy among the neutrino families. In this case, if neutrino masses m v S3 30 
eV, they could still contribute significantly to the non-baryonic dark matter in the 
Universe. 

These neutrinos left over from the Big Bang were in thermal equilibrium in the 
early Universe and decoupled when they were moving with relativistic velocities. 
They fill the Universe in enormous quantities and their current number density is 
similar to the one of photons (by entropy conservation in the adiabatic expansion of 
the Universe). In particular, n v = A- n 1 . Moreover, the number density of photons 
is very accurately obtained from the cosmic microwave background measurements. 
The present temperature T S3 2.725 K implies n 7 S3 410.5 cm' 3 . Thus one can 
compute the neutrino mass density p v = mtot n v , where m to t is basically the total 
mass due to all flavours of neutrino. Hence, 

VL V ss . (4) 

30 eV y ' 

Clearly, neutrinos with m v < 1 eV cannot solve the dark matter problem, but a 

neutrino with m v S3 30 eV would give £7^ ss 1 solving it. Unfortunately, due to the 

small energies involved, detection of these cosmological neutrinos in the laboratory 

is not possible. 

However, even if m v ss 30 eV, there is now significant evidence against neutrinos 
as the bulk of the dark matter. Neutrinos belong to the so-called 'hot' dark matter 
because they were moving with relativistic velocities at the time the galaxies started 
to form. But hot dark matter cannot reproduce correctly the observed structure 



February 1, 2008 22:9 WSPC/INSTRUCTION FILE 

darkijmpa21+refs+refs+refs 



Dark Matter Detection in the Light of Recent Experimental Results 9 

in the Universe. A Universe dominated by neutrinos would form large structures 
first, and the small structures later by fragmentation of the larger objects. Such a 
Universe would produced a 'top-down' cosmology, in which the galaxies form last 
and quite recently. This time scale seems incompatible with our present ideas of 
galaxy evolution. 

This lead to fade away the initial enthusiasm for a neutrino-dominated Uni- 
verse. Hence, many cosmologists now favour an alternative model, one in which the 
particles dominating the Universe are 'cold' (non-relativistic) rather than hot. This 
is the case of the axions and neutralinos which we will study below. 



3.1.2. Axions 

Axions are spin particles with zero charge associated with the spontaneous break- 
ing of the global U(l) Peccei-Quinn symmetry, which was introduced to dynami- 
cally solve the strong CP problem ""4 Let us recall that this important problem 
arises because QCD includes in its Lagrangian the term j^tr{F liV F' lv ) violating 
CP, and therefore there are important experimental bounds against it. In particu- 
lar, the stringent upper limit on neutron dipole electric moment implies the bound 
6 < 10~ 9 . This is the problem. Why is this value so small, when a strong interaction 
parameter would be expected to be 0(1)? 

Although axions are massless at the classical level they pick up a small mass by 
non-perturbative effects. The mass of the axion, m a , and its couplings to ordinary 
matter, g a , are proportional to l/f a , where f a is the (dimensionful) axion decay 
constant which is related to the scale of the symmetry breaking. In particular, 
the coupling of an axion with two fermions of mass my, is given by g a ~ m/// a . 
Likewise, m a ~ A 2 CD / '/„, i.e. 

„^ 10 -eVxi2^SX. (5, 
J a 

A lower bound on f a can be obtained from the requirement that axion emission 
does not over-cool stars. The supernova SN1987 put the strongest bound, f a > 10 9 
GeV. On the other hand, since coherent oscillations of the axion around the min- 
imum of its potential may give an important contribution to the energy density 
of the Universe, the requirement O < 1 puts a lower bound on the axion mass 
implying f a < 10 12 GeV. The combination of both constraints, astrophysical and 
cosmological, give rise to the following window for the value of the axion constant 

10 9 GeV < /„ < 10 12 GeV . (6) 

The lower bound implies an extremely small coupling of the axion to ordinary 
matter and therefore a very large lifetime, larger than the age of the Universe by 
many orde rs of magnitude. As a consequence, the axion is a candidate for dark 
matter ESI. Axions would have been produced copiously in the Big Bang, they were 
never in thermal equilibrium and are always nonrelativistic (i.e. they are cold dark 
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matter). In addition the upper bound implies that m a ~ 10 -5 eV if the axion is to 
be a significant component of the dark matter. 

Since the axion can couple to two photons via fermion vacuum loops, a pos- 
sible way to detect it is through conversion to photon in external magnetic field. 
Unfortunately, due to the small couplings to matter discussed above, we will not 
be able to produce axions in the laboratory. On the other hand, relic axions are 
very abundant (as we will discuss in Section 4, the density of dark matter particles 
around the Earth is about 0.3 GeV cm -3 , since m a ~ 10 -5 eV there will be about 
10 13 axions per cubic centimeter) and several experiments are trying already to 
detect axions or are in project. For example, an experiment at Lawrence Livermore 
National Laboratory (California, US) has reported in 1998 its first results excluding 
a particular kind of axion of mass 2.9 x 10 -6 eV to 3.3 x 10 -6 eV as the dark matter 
in the halo of our galaxy. 

Let us finally mention that links to web pages of axion detection experiments 



can be found in the web page http://cdms.physics.ucsb.edu/others/othcrs.html 



3.1.3. WIMPs 

As mentioned in the Introduction, WIMPs are very interesting candidates for dark 
matter in the Universe. They were in thermal equilibrium with the standard model 
particles in the early Universe, and decoupled when they were non-relativistic. The 
process was the following. When the temperature T of the Universe was larger than 
the mass of the WIMP, the number density of WIMPs and photons was roughly 
the same, n cx T 3 , and the WIMP was annihilating with its own antiparticle into 
lighter particles and vice versa. However, after the temperature dropped below the 
mass of the WIMP, m, its number density dropped exponentially, n cx e -m / T , 
because only a small fraction of the light particles mentioned above had sufficient 
kinetic energy to create WIMPs. As a consequence, the WIMP annihilation rate 
r = ((T ann v)n dropped below the expansion rate of the Universe, T < H. At this 
point WIMPs came away, they could not annihilate, and their density is the same 
since then (freeze-out typically occurs at Tp as m/20). This can be obtained using 
the Boltzmann equation, which describes the time evolution of the number density 
n(t) of WIMPs 

^+3Hn = -(a ann v)[(n) 2 -(n^) 2 ] , (7) 

where H is the Hubble expansion rate, cr ann is the total cross section for annihilation 
of a pair of WIMPs into standard model particles, v is the relative velocity between 
the two WIMPs, (...) denotes thermal averaging, and n eq is the number density of 
WIMPs in thermal equilibrium. One can discuss qualitatively the solution using 
the freeze-out condition T = ((J a nnv) F n = H . Then, the current WIMP O/i 2 = 
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Fig. 2. Feynman diagrams contributing to early-Universe neutralino (x?) annihilation into 
fermions through neutral Higgses (H = H,h,A) and squarks and sleptons (/). 



(p/p c )h 2 , turns out to be 

ntf * 3 x 10 ' 27 cm3 8-1 , (8) 

(farm V) 

where the number in the numerator is obtained using the value of the temperature 
of the cosmic background radiation, the Newton's constant, etc. As expected from 
the above discussion about the early Universe, the relic WIMP density decreases 
with increasing annihilation cross section a . 

Now we can understand easily why WIMPs are so good candidates for dark 
matter. If a new particle with weak interactions exists in Nature, its cross section 
will be a ~ a 2 /m 2 oak , where a ~ 0(1O~ 2 ) is the weak coupling and m wcak ~ 0(100 
GeV) is a mass of the order of the W gauge boson. Thus one obtains a 10~ 9 
GcV~ 2 w lpb (recall that in natural units 1 GeV" 2 = 0.389 x 10~ 27 cm 2 = 
0.389 x 10 9 pb). Since at the freeze-out temperature Tp the velocity v is a significant 
fraction of the speed of light (v 2 sa c 2 /20), one obtains < a an n v >» 10~ 26 cm 3 
s . Remarkably, this number is close to the value that we need in Eq. JSJ in order 
to obtain the observed density of the Universe. This is a possible hint that new 
physics at the weak scale provides us with a reliable solution to the dark matter 
problem, and also a qualitative improvement with respect to the axion dark matter 
case, where a small mass for the axion about 10~ 5 eV has to be postulated. 

SUSY is a theory that introduces new physics precisely at the weak scale, and 
that, as discussed in the Introduction, predict a new particle, the neutralino, which 
could be stable. These are the reasons to consider the neutralino as a very seri- 
ous candidate for the sought-after dark matter. Concerning the annihilation cross 
section contributing to the density of the Universe in Eq. (JHJ) , there are numer- 
ous final states into which the neutralino can annihilate. The most important of 
these are the two body final states which occ ur at the tree level. Specially these 
are fermion-antifermion pairs 1 W l 1 Q I I (// where / are the standard model 

a Let us remark that the theoretical computation of the relic density depends on assumptions 
about the evolution of the early Universe, and therefore cosmological scenarios different from the 
thermal production of neutralinos discussed here would give rise to different results I2Z| 
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quark s and leptons), as those shown in Fig. [2J Others are wea k gauge bosons pairs 
l lCTWUfl W + W~, Z°Z°, and those containing Higgs bosons E ETO I lOfllU ll suc h 
as W+H-, W~H+, Z°A, Z°H, Z°h, H + H~ and all six combinations of A,h,H. 
Different subtleties of the analyses has been discussed in Refs. llOU 11051 The anni- 
hilation cross section is of the form 

where Aann is the amplitude which depends on the dynamics of the collision and 
N ann is the number of annihilation channels. Generically A an n ~ 1/A4 2 , where 
M is the mass of the particles mediating the interaction. Thus (7 ann ~ l/m?~ 0j 
and therefore in order to satisfy the upper bound in Eq. J2J, moderate values of 
the LSP mass are necessary. Several exceptions to this rule will be discussed in 
Subsection 5.2.1. Let us mention here one of them which is particularly interesting. 

In principle to use the above discussed neutralino annihilation cross section 
is sufficient, since the effects of heavier particles are suppressed by the Boltz- 
mann factor. However, the next to lightest supersymmetric particle (NLSP) may 
lie near in mass to the LSP so that both particles freeze out of equilibrium 
at approximately the same temperature. Thus the NLSP should be included in 
principle in the reaction network, since coannihilatio n ch annels NLSP-LSP (and 
also channels NLSP-NLSP) might be now relevant EHI l n fact, this is so only 
when (m NLSP — m^o)/m^o < 0.2, since the NLSP number density is suppressed 

by e ^ mNLSP m x?) Tp re i a tive to the neutralino number density, where recall that 
Tp w m^o/20. We will see in Section 5 that special regions of the parameter space of 
the MSSM fulfil this condition, and therefore including coannihilations is important. 
Coannihilation channels in SUSY, in particular Yi-stau, Y?-charginq and _Y? -stop, 
have been exhaustively studied in the literature llE 6 | 107 | l08 | 109 | ilO | lll | 112 r 

Taking into account all the above comments concerning SUSY scenarios, we 
will see in detail in Section 5 that significant regions of the parameter space of the 
MSSM produce values of the cross section in the interesting range mentioned below 
Eq. ©. 

On the other hand, since neutralinos, or WIMPs in general, interact with or- 
dinary matter with roughly weak strength, their presence in galactic scales, and 
in particular in our galaxy, raises the hope of detecting relic WIMPs directly in a 
detector by observing their scattering on target nuclei through nuclear recoils. This 
will be the subject of the next section, but before analyzing it let us briefly discuss 
here other candidates for dark matter. 



3.2. Other candidates 

The first proposed candidate for dark matter in the co ntext of WIMPs was in fact a 
heavy (Dirac or Majorana) fourth-generation neutrino ^^'^ with the usual standard 
model couplings. It is not easy to construct theoretical models with such a stable 
heavy particle, but in any case LEP limits from the Z° width imply to„ > m z o/2, 
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and this is not compatible with a correct relic abundance. The relic abundance can 
be studied as in Subsection 3.1.3 with the annihilation of the heavy neutrinos into 
light fermions through the exchange of Z . The final result i s that ne utrinos with 
the a bove l ower bound give rise to a too small relic abundance ^E^Sl j n addition, 
direct and indirect Illfil detection experiments rule out m u < 1 TeV. Attempts 
to consider still massive fourth-generation neutrinos as dark matter candidates, can 
be found in Ref. ITTTI 

In addition to the neutralino there are other potential SUSY candidates for dark 
matt er. In principle one of them mig ht be the sneutrino^^ of the MSSM. One 
finds 111*3 that the sneutrino relic density is in the region 0.1 < fl^h 2 < 1.0 for 550 
GeV <ma< 2300 GeV. This is consistent with the LEP limits on Z° -> invisible 
neutral particles suggesting mp > m z o /2, as in the case of the heavy neutrino above. 
However, sneutrino-nucleus interaction is similar to t he he avy neutrino-nucleus in- 
teraction, and therefore direct detection experiments E5J impose similar limits on 
the sneutrino mass as above. Such a heavy sneutrino cannot be the LSP in SUSY 
models. 

The supersymmetric partner of t he gr aviton, the gravitino, has also been pro- 
posed as a candidate for dark matter 11*121. In the absence of inflation, it could give 
rise to the correct relic abundance if its mass is of order keV. This is unlikely in 
specific theoretical models. For example, in gravity mediated SUSY breaking the 
masses of the superpartners are of order the gravitino mass, and therefore this 
should be of order 1 TeV. Gravitinos as dark matter would be undetectable since 
their interactions with ordinary matter are extremely weak. 

Another candidate for dark matter that can arise in SUSY is the cmnoEIl 
Models combining SUSY and the Peccei-Quinn solution to the strong CP problem 
necessarily contain this particle. Unlike the neutralino and gravitino, the axino 
mass is generically not of order the SUSY-breaking scale and can be much smaller. 
For a recent work with masses in the range of tens of MeV to several GeV as 
corresponding to cold axino relics see e.g. Ref. 11221 

The possibility that dark matter coul d be made of additional scalars in the stan- 
dard model has also been analyzed E23I Scalar ca ndid ates have also been studied 
in a prototypical theory space 'Little Higgs Model' Ell a nd in theories with N = 2 
extended supersymmetry and/or extra space dimensions Ell 

More (exotic) candidates can be found in the literature. Let us mention 
some of them . Strongly Interacting Massive Particles, SIMPs, have been proposed 
| 12fi | 127 |12g| p Qr example, bound state of ordinary quarks and gluons with a heavy 
stable quark, scalar quark, gluino, or colored technibaryon. 

Mirror matter is predicted to exist if Nature exhibits a n exa ct unbroken mirror 
symmetry, and then could be a candidate for dark matter LL2*2I 

D ark m atter made of CHArged Massive Particles, CHAMPs, has also been pro- 
posed ESS Apparently, if their masses are larger than 20 TeV these particles will 
very rarely be trapped in ordinary matter, and are safe from bounds on exotic 
heavy isotopes. 
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A potential conflict between collisionless dark matter scenarios and observations 
motivated the authors of Ref. I Kill to propose Self Interacting Dark Matter, SIDM, 
i.e. dark matter with a large scattering cross section but negligible annihilation. 
A concrete realization of this idea, using a scalar gauge singlet, can be found in 
Ref. 11321 On the other hand, the authors of Ref. 11331 solved the problem using a 
mechanism of non-thermal production of WIMPs. 

Instead of using WIMPs with typical masses of order a hundred GeV, the au- 
thors of Ref. 11341 studied scenarios with nonthermal WIMPs in the range 10 12 to 
10 16 GeV. They called these objects WIMPZILLAS. 

The lightest Kaluza-Klein particle, LKP, in model s wit h TeV extra dimensions 
has been studied as a viable dark matter candidate E3. It is actually a typical 
WIMP (the most studied possibility is a Kaluza-Klein 'photon'), with a mass in 
the range 400-1200 GeV. 

Superweakly-interacting massive particles, superWIMPs, were also proposed 
in Ref. 11361 They naturally inherit the desired relic density from late decays of 
metastable WIMPs. Well-motivated examples are weak-scale gravitinos in super- 
gravity and Kaluza-Klein gravitons from extra dimensions. 

In the context of the heterotic string Eg x Eg, where there is a natural hidden 
sector associated with the second Eg, it was mentioned in Ref. 11371 that this sector 
which only interacts with ordinary matter through gravitational interactions could 
be a shadow world. The astrophysical and cosmological implications of this comment 
were analyzed in detail in Ref. 138 

Still, in the context of strings, other possibilities were also analyzed. For exam- 
ple, in Ref. 11391 the authors proposed Cryptons, which are stable bound states of 
matter in the hidden sector. In Ref. 140 glucballs of the hidden-sector non-Abclian 
gauge group were proposed as candidates for SIDM. Brane-world dark matter was 
studied in Ref. 141, 

Let us finally mention that the list of dark matter candidates discussed here is 
large but by no means complete. The reason being that different candidates have 
been proposed during many years, and therefore it is almost impossible not to forget 
some of them. In fact, this list is still increasing (see e.g. the works in Ref. I142|) . 

4. Dark matter detection 

Given the discussion in the previous section, one can say that there are good particle 
candidates for dark matter. As a matter of fact, we saw that WIMPs and axions 
are particularly interesting. Since the former can be left over from the Big Bang in 
sufficient number to account for the relic matter density in a natural way, we will 
concentrate on them. Thus we will review in this section the current and projected 
experiments for detecting WIMPs (e.g. the detection of axions was briefly discussed 
in Section 3.1.2). In particular we will analyze the two possible ways of detecting 
them: direct and indirect detection. Of course, the most clear one is the direct 
detection by observing their scattering on the material in a detector, and we will 
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Fig. 3. Elastic scattering of a dark matter particle with an atomic nucleus in a detector. 

study it in detail. Finally, we will discuss briefly the indirect detection by looking 
for evidence of their annihilation. 

4.1. Direct detection 

As discussed in the Introduction, if neutralinos, or WIMPs in general, are the bulk 
of the dark matter, they will form not only a background density in the Universe, 
but also will cluster gravitationally with ordinary stars in the galactic halos. In 
particular they will be present in our own galaxy, the Milky Way. This raises the 
hope of detecting relic WIMPs directly, by experiments on the Earth. In particular, 
through scattering with the material in a detector. In fact general studies of the 
possibility of dark matter detection began around 1982. Since the detection will be 
on the Earth we need to know the properties of our galaxy in order to be sure that 
such a detection is feasible b . 

As a matter of fact, rotation curves are much better known for external galaxies 
than for ours, due to the position of the Earth inside the galaxy. In any case, 
analyses have been carried out with t he co nclusion that indeed the Milky Way 
contains large amounts of dark matter EH. Besides, some observational evidence 
seems to point at a roughly spherical distribution of dark matter in the galaxy. At 
the position of the Sun, around 8.5 kpc away from the galactic center, the mean 
density of elementary particles trapped in the gravitational potential well of the 
galaxy is expected to be p ~ 5 x 10" 25 gr cm" 3 ~ 0.3 GeV cm" 3 . For WIMPs 
with masses about 100 GeV this means a number density tiq w 3 x 10" 3 cm" 3 . In 
addition, their velocity will be similar to the one of the Sun since they move in the 
same gravitational potential well, vq « 220 km s , implying a flux of dark matter 
particles Jo = uq vq w 10 5 cm" 2 s" 1 incident on the Earth c . Although this number 

b Direct detection of extragalactic WIMPs has also been analyzed recentlyESl Although of much 
lower flux than typical galactic halo WIMPs, it seems that they have a number of interesting 
features for their detectability. 

c Clumps of dark matter crossing the Earth with a density larger by a factor 10 s , and with a 
periodicity of 30-100 Myrs, have been proposed to explain double mass extinctions of paleontology. 
First, the direct passage of a dark matter clump may lead to a preliminary extinction 

145 by 
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Fig. 4. Feynman diagrams contributing to neutralino-nucleon cross section through squark (q) 
exchange and CP-even light (h) and heavy (H) neutral Higgs exchange. 

is apparently large, the fact that WIMPs interact weakly with matter makes their 
detection very difficult. Most of t hem will pass through matter without prevention. 
In any case, as suggested in 1985 -33 direct experimental detection of WIMPs is in 
principle possible. The two types of direct detection experiments will be discussed 
below. 

(i) Elastic scattering 

The detection of WIMPs through elastic scattering with nuclei in a detector is 
shown schematically in Fig. [21 As we can see the nucleus recoils as a whole. A very 
rough estimate of the rate R in a detector is the following. A particle with a mass 
of order 100 GeV and electroweak interactions will have a cross section er « 1 pb, 
as discussed in Subsection 3.1.3. Thus for a material with nuclei composed of about 
100 nucleons, i.e. M N ~ 100 GeV = 177 x 10~ 27 kg, one obtains R ~ J a/M N w 10 
events kg -1 yr . This means that every day a few WIMPs, the precise number 
depending on the number of kilograms of material, will hit an atomic nucleus in a 
detector. 

Of course the above computation is just an estimate and one should take into 
account in the exact computation the interactions of WIMPs with quarks and 
gluons, the translation of these into interactions with nucleons, and finally the 
translation of the latter into interactions with nuclei. In the case of ncutralinos 
as WIMPs, diagrams contributing to neutralino-quark cross section are shown in 
Fig. 0] The relevant (scalar) xj-nucleus cross section is of the form 

a scat ~ M?\A sca t\ 2 , (10) 

where M r — M^m^o /(-Mjv + rn^o) is the reduced mass with Mn the mass of the 
nucleus, and A sca t is the amplitude which depends on the dynamics of the collision. 

causing lethal carcinogenesis in organisms ESI Later, the accumulation of dark matter in the 
center of the Earth and the subsequent annihilation would produce a large amount of heat with 
the consequent ejection o f sup erplumes . followed by massive volcanism leading to the second burst 
of more severe extinction! ' I For another exotic mechanism solving mass extinctions, and related 
to dark matter, see Ref. 11481 where mirror matter is used. 



February 1, 2008 22:9 WSPC/INSTRUCTION FILE 

darkijmpa21+refs+refs+refs 



Dark Matter Detection in the Light of Recent Experimental Results 17 

In particular, the quarks masses m qi the hadronic matrix elements ff? , the proton 
mass m p , and the masses of the particles mediating the interaction, such as rriq, 
nriHi m-h, enter in A SC at- We will check in Section 5 that significant regions of the 
parameter space of the MSSM produce values of the neutralino-nucleus cross section 
a sca t — 1 pb, and therefore giving rise to a reasonable number of events. As a matter 
of fact, in the experimental results that one finds in the literature the authors prefer 
to give the WIMP-nucleon cross section. As will be discussed in Section 5, this is 
about eight orders of magnitude smaller than the WIMP-nucleus cross section, and 
therefore a typical value is w 10 -8 pb. In the next Sections, when talking about 
scattering cross section, we will always consider this one. 

Let us finally remark that the diagrams for neutralino annihilation (see Fig. [5J 
are related to these by crossing symmetry. Thus, provided that the main annihi- 
lation channel is into fermions, the amplitudes of annihilation and scattering with 
nucleons are similar, and this leads for the amplitudes with the nucleus to the rela- 
tion |*4 S cat| 2 — c 2 A 2 \A an n\ 2 , where A is the atomic weight and c 2 is a constant (we 
can deduce from Section 5 that c ~ fjl m p /m s = 0(1)). From Eqs. JSJ and 110JI 
it is obvious that o~ sca tlo~ ann ~ const. However, if the neutralinos are heavy they 
have other annihilation channels, such as Higgs bosons or vector boson pairs, and 
therefore the crossing argument does not apply. 

(ii) Inelastic scattering 

There are two possible ways of detecting WIMPs directly through inelast ic sca tter- 
ing. These can be with nuclei in a detector ESI or with orbital electrons ESI The 
former produces an excited nuclear state, and has the double signature of recoil 
energy followed ~ 1 ns later by the emission of a decay photon. However, the event 
rates seem to be suppressed and problems of natural radioactivity and expense 
disfavour the materials that might be used for the detector d . A similar conclu- 
sion concerning the event rates was obtained for the case of scattering of WIMPs 
from orbital electrons. These could leave an excited state emitting a photon by the 
deexcitation. 

4.1.1. Experiments around the world 

More than 20 experiments for the direct detection of dark matter are running 
or in preparation around the world. Given the above results, all of them use the 
elastic-scattering technique. For example, Germanium is a very pure material and 
has been used for many years for detecting dark matter in this way. In this type of 
experiments, in order to detect the nuclear recoil energy, they measure the ionization 
produced by collision with electrons. In f act, 76 Ge ionization detectors has been 
applied to WIMP searches since 1987ESI In 2000 the situation was the following. 



d For a recent analysis with a different conclusion, see Ref . 11511 
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Fig. 5. Earth's motion around the Sun. 

The best combination of dat a from these experiments, to gethe r with the last data 
from the Heidelberg-Moscow and IGEX experiments l ^^ l located at the Gran 
Sasso (L'Aquila, Italy) and Canfranc (Huesca, Spain) Underground Laboratories, 
respectively, were able to exclude a WIMP-nucleon cross section larger than about 
10 -5 pb for masses of WIMPs ~ 100 GeV, due to the negative search result. 
Although this was a very interesting bound, it was still well above the expected 
value ~ 10~ 8 pb. 

Let us remark that it is convenient to carry the experiments out in the deep 
underground. For a slow moving (~ 300 km s -1 ) and heavy (~ 100 — 1000 GeV) 
particle forming the dark matter halo, the kinetic energy is very small, around 
100 keV, and in fact the largest recoil energy transfer to a nucleus in the detector 
will only be a few keV. Since cosmic rays with energies ~ keV-MeV bombard the 
surface of the Earth, the experiments must have an extremely good background 
discrimination. In particular, neutrons coming from collisions between cosmic-ray 
muons and nuclei produce nuclear recoils similar to those expected from WIMPs at 
a rate ~ 10 3 events kg -1 day -1 . Thus detectors located in the deep underground, 
reduce the background by orders of magnitude in comparison with the sea level 
intensity. 

In fact, this is still not enough since the detector has to be protected also against 
the natural radioactivity from the surroundings (e.g. the rocks) and the materials of 
the detector itself. This produces again neutrons but also X rays, gamma rays and 
beta rays giving rise to electron recoils. The latter may be a problem for detectors 
based only on ionization or scintillation light since nuclear recoils with energies of 
a few keV are much less efficient in ionizing or giving light than electrons of the 
same energy. Various protections aim to reduce these backgrounds. In particular, 
low radioactive materials, such as e.g. high-purity copper or aged lead, are used for 
the shielding. In addition, high-purity materials for the detector are also used. 

Summarizing, with this type of experiments the WIMP nuclear recoil signal will 
appear as an excess of recoil events above the expected background rate. However, 
it would be very interesting to also look for some additional feature of the WIMP 
signal that positively identifies it as galactic in origin. In this sense a different 
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method for discriminating a dark matter signal from background is the so-called 
annual modulation signature l l^ h As it is shown schematically in Fig. El as the 
Sun orbits the galaxy with velocity vq rs 220 km s _1 , the Earth orbits the Sun with 
velocity w 30 km s _1 and with the orbital plane inclined at an angle of 60° to the 
galactic plane. Thus e.g. in June the Earth's rotation velocity adds to the Sun's 
velocity through the halo with a maximum around June 2, whereas in December 
the two velocities are in opposite directions. When this is taken into account the 
Earth velocity is given by 

2n(t — t m 
1 year 

where t m = June 2± 1.3 days. This fluctuation produces a rate variation between 
the two extreme conditions. The variation is so small « 7% that the experiment can 
only work if large number of events are found, implying that large mass apparata 
are necessary. 

The DArk MAtter (DAMA) experiment investigates the annual modulation 
of this signature rather than discriminating signal events against background. It 
consists of about 100 kg of material in a small room at the Gran Sasso National 
Laboratory located beside the Gran Sasso Tunnel on the highway connecting Ter- 
amo to Rome (see Fig. |§] for a similar experiment DAM A/LIBRA). The maximum 
thickness of the rock overburden is 1400 m. In the experiment they use nine 9.70 
kg Nal crystal scintillators which measure the ionization produced by the nuclear 
recoil through the emission of photons. Remarkably, they found that their detec- 
tors flashed more times in June than in December. The data collected El over four 
yearly cycles, DAMA/NaI-1,2,3,4, until the second half of August 1999, strongly 
favour the presence of a yearly modulation of the rate of the events. Moreover, very 
recently, more data ha v ; been reported E I confirming this result. In particular the 
DAMA/NaI-5 data have been collected from August 1999 to end of July 2000. Af- 
terwards, the DAMA/NaI-6 data have been collected from November 2000 to end 
of July 2001, while the DAMA/NaI-7 data have been collected from August 2001 to 
July 2002. The analysis of the data of the seven annual cycles offers an immediate 
evidence of the presence of an annual modulation of the rate of the events in the 
lowest energy region as shown in Fig. (from Ref. 0, where the time behaviours 
of the (2-4), (2-5) and (2-6) keV residual rates are depicted. This signal is com- 
patible 1^1 with WIMP masses up to 100 GeV and WIMP-nucleon cross sections in 
the i nterv al 10 -6 — 10 -5 pb, as shown with a dark shaded area in Fig. [TJ where the 
data HSi3 from DAMA/NaI-0 have also been taken into account. The lower bound 
of 30 GeV in the figure is purely theoretical, and inspired by the lower bound on 
the lightest neutralino of SUSY models, as derived from the LEP data using GUT 
conditions. 

c Another novel possibility to search for characteristic signatures has recently been proposed, 
namely by obse rving the electrons which follow the ionization of the atom during the WIMP- 
nucleus collision HSU 
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Fig. 6. Model independent residual rate for events in the (2—4), (2—5) and (2-6) keV energy 
intervals as a function of the time elapsed since January 1-st of the first year of data taking. 
The experimental points present the errors as vertical bars and the associated time bin width as 
horizontal bars. The superimposed curves represent the cosinusoidal functions behaviours expected 
for a WIMP signal with a period equal to 1 year and phase at 2 nd June. 



It is worth remarking that this result has been obtained assuming the simple 
isothermal sphere halo model with dark-matter density po = 0.3 GeV cm -3 and 
WIMP velocity v = 220 km s" 1 (to be precise one assumes a Maxwcll-Boltzmann 
local velocity distribution f(v) producing a velocity dispersion v =< 

v 2 > 1 ' 2 = (3/2) 1 / 2 vo ~ 270 km s _1 ). However, when uncertainties on the halo 
model are taking into account, the signal is consistent with a larger region of the 
parameter space. In particular, in Refs. E] and 11571 modifications in the velocity 
distribution function for different galactic halo models were considered, using in 
addition the allowed ranges for vq and po in each model (see Ref. 11581 for other 
studies of uncertainties on the halo model). The final result of the analyses is shown 
in Fig.0with a light shaded area. One sees that the signal is compatible with larger 
values of the parameters (for a different opinion see Ref. I159| ). i.e. WIMP masses 
up to 270 GeV, and WIMP-nucleon cross sections in the interval 10~ 7 — 6 x 10~ 5 
pb. This result corresponds at a rate of about 1 event per kg per day. In fact, as 
discussed in Ref. 11571 (see also Ref. |SJ) , when co-rotation of the galactic halo is also 
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Fig. 7. Areas allowed by the different experiments for the direct detection of dark matter in 
the parameter space ((JwiMP-nucicon, TOwimp), where uwiMP-nucieon is the WIMP-nucleon cross 
section and m W i M p is the WfMP mass. The light shaded (yellow) area is allowed by the DAMA 
experiment, including uncertainties on the halo model. The dark shaded (brown) area inside the 
previous one corresponds to the simple case of an isothermal sphere halo model. This case is 
also the one analyzed by the other experiments: The (lower) areas bounded by solid, dot-dashed, 
dashed and (again) solid lines are allowed by IGEX, CDMS, ZEPLIN I and EDELWEISS current 
experimental limits. The (upper) areas bounded by solid, dot-dashed and dashed lines will be 
analyzed by the projected GEDEON, CDMS Soudan and GENIUS experiments. 



considered, the mass range extends further to 500 — 900 GeV, for cross sections in 
the interval fewxl0~ 6 — 2 x 10~ 5 pb. 

Although the DAMA group is confident about the data, since the y claim to have 
ruled out systematic effects which could fake the signature H50HBI I as e g temper- 
ature changes, it is worth remarking that the above values for the cross section are 
generically above the expected weak-interaction value, and therefore they are not 
easy to obtain in SUSY models with neutralino dark matter, as we will discuss in 
Section 5. But in fact, the DAMA result is controversial, mainly because the nega- 
tive search result obtained by other recent experiments. The first of these was the 
Cryogenic Dark Matter Search (CDMS) experiment ® in the US. This was located 
just 10 metres below ground at Stanford University in California, and therefore it 
must discriminate WIMPs signals against interactions of background particles. Two 
detection techniques are used for this discrimination, both the ionization and the 
temperature rise produced during a recoil are measured. The latter can be observed 
since the recoiling nucleus is stopped within 10~ 7 — 10~ 6 cm (~ 10~ 14 s) releas- 
ing a spherical wave of phonons traveling at ~ 5 x 10 5 cm s _1 , and subsequently 
converted to a thermal distribution. These two techniques allow to discriminate 
electron recoils caused by interactions of background particles from WIMP-induced 
nuclear recoils. The ratio of deposited energies heat/ionization would be ~ 2 — 3 
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for the former and larger than 10 for the latter. However, although neutrons are 
moderated by a 25-cm thickness of polyethylene between the outer lead shield and 
cryostat, an unknown fraction of them still survives. Two data sets were used in 
this analysis: one consisting of 33 live days taken with a 100-g Si detector between 
April and July 1998, and another consisting of 96 live days taken with three 165-g 
Ge detectors between November 1998 and September 1999. Although four nuclear 
recoils are observed in the Si data set, they cannot be due to WIMPs, they are due 
in fact to the unvetoed neutron background. On the other hand, in the Ge detec- 
tor thirteen nuclear recoils are observed in the 10.6 kg per day exposure between 
10 and 100 keV, which is a similar rate to that expected from the WIMP signal 
claimed by DAMA. However, the CDMS group concludes that they are also due to 
neutrons. These data exclude (see also Ref. 11621 for more recent data) much of the 
region allowed by the DAMA annual modulation signal f , as shown in Fig. s . 

In addition, a small part of the region excluded by CDMS has also been excluded 
by IGEX 1— , as shown in Fig. [71 To obtain this result, 40 days of data from the 
IGEX detector were analyzed. Let us recall that this experiment is located at the 
Canfranc Tunnel Astroparticle Laboratory in the international railway tunnel of 
Somport at Canfranc (Spanish pyrenees), under a 860 m rock overburden, and 
that consists of 2 kg Germanium ionization detector. A region similar to the one 
excluded by IGEX h as al so been excluded by the Heidelberg Dark Matter Search 
(HDMS) experiment Ell This operates two ionization Ge detectors in a unique 
configuration, and was installed at Gran Sasso in August 2000. The data used for 
the analysis were taken from February 2001 to September 2001. 

But more disturbing are the recent results from EDELWEISS and ZEPLIN I 
collaborations excluding even larger regions than CDMS, as shown also in Fig. [TJ 
The EDELWEISS experimental is located at the Frejus Underground Laboratory in 
the Frejus Tunnel under the French-Italians Alps, under a 1780 m rock overburden. 
As CDMS, this experiment also uses a heat-and-ionization cryogenic Ge detector, 
although they differ by their mass, geometry and electrode implantation scheme. 
First, this collaboration was using a 320 g detector, 70 mm in diameter and 20 mm 
in height, but this was not sufficient to extend the sensitivity to the values obtained 
by DAMA. Following this result, three new detectors with the same characteristics 
as the first one were put in operation simultaneously. The data consists of all physics 
runs recorded over a period from February to May 2002. 

On the other hand, the ZEPLIN project^ at the Boulby salt mine (Yorkshire, 
UK) consists of a series of liquid Xenon detectors operating some 1100 m under- 

'For some attempts to show that both experiments, DAMA and CDMS, might not be in conflict, 
see Ref. EH 

g As mentioned in the caption of Fig. [7] unlike DAMA for CDMS and for the other experiments 
analyses taking into account the uncertainties in the galactic halo are not shown and we only 
see the effect of the standard halo model on their results. Including those uncertainties, the light 
shaded area favoured by DAMA and not excluded by the null searches would be in principle 
smaller than the one shown here (for an analysis of this see Ref. 11(541 . 
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Fig. 8. Left picture: during the LIBRA detectors installation in HP Nitrogen atmosphere. Right 
picture: view at end of the detectors installation. All the used materials have been deeply selected 
for radiopurity (see for example the cables with teflon envelop). 



ground, where the nuclear recoil produces both an ionization and a scintillation 
signal, potentially capable of providing much greater signal discrimination than 
Nal. The ZEPLIN I detector has been running for approximately one year with a 
mass of about 4 kg. The collaboration is concentrating on new developments for 
future ZEPLIN detectors. 

4.1.2. Future dark matter searches 

Owing to this controversy between DAMA and the other experiments, one cannot 
be sure whether or not the first direct evidence for the existence of dark matter has 
already been observed. Fortunately, the complete DAMA region will be tested by 
current dark matter detectors. This is for example the case of the IGEX experi ment 
mentioned above, with a moderate improvement of the detector performances^^. 
In addition, a new expe rimen tal project, GErmanium DEtectors in ONe cryostat 
(GEDEON), is planned 023. i t w m use the technology developed for the IGEX 
experiment and it would consist of a set of ~ 1 kg Germanium crystals, of a total 
mass of about 28 kg, placed together in a compact structure inside one only cryostat. 
GEDEON would be massive enough to search for the annual modulation effect 
and explore positively a WIMP-nucleon cross section a > 3 x 10~ 8 pb (see Fig. EJ). 
Likewise an improved HDMS e xper iment should be able to test the DAMA evidence 
region within about two years ' 

DAMA and CDMS collaborations also plan to expand their experiments. The 
DAMA collaboration dismounted the 100 kg Nal set-up and installed the new LI- 
BRA set-up. The latter consisting of about 250 kg of Nal made of 25 detectors, 9.70 
kg each one. This will make the experiment more sensitive to the annual modula- 
tion signal. Some pictures taken during the installation can be seen in Fig. [S] (from 
Ref. UnTJ). This was completed at the end of 2002. On the other hand, the CDMS 
collaboration is moving its detector to the abandoned Soudan mine in Minnesota 
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^^n^^!^^SS^^^^^^^^^^H^ Liquid Nitrogen Isolation Steel-vessel 

(a) ^Boron Polyethelene (b) 



Fig. 9. (a): Conceptual design of the Genius TF. Up to 14 detectors will be housed in the inner 
detector chamber, filled with liquid nitrogen. As a first shield 5 cm of zone refined Germanium, or 
extremely low-level copper will be used. Behind the 20 cm of polystyrene isolation another 35 cm 
of low level lead and a 15 cm borated polyethylene shield will complete the setup, (b): GENIUS - 
100 kg of Ge detectors are suspended in a large liquid nitrogen tank. 

(approximately 700 metres below ground), increasing also the mass of its Ge/Si 
targets to about 10 kg by 2006. This experiment, CDMS Soudan, will be able to 
test a WIMP-nucleon cross section a > 10~ 8 pb (see Fig. [TJ) . 

But this is not the end of the story, since a new generation of very sensitive 
experiments have been proposed all over the world. For instance, only in the Gran 
Sasso laboratory there will be five experiments searching for WIMPs. Apart from 
the two already discussed DAMA and HDMS, there are three other experiments 
in prospect, CRESST, CUORE and GENIUS. For example, the Cryogenic Rare 
Event Search using Superconducting Thermometers (CRESST) experiment EBHI 
measures simultaneously phonons and scintillation light distinguishing the nuclear 
recoils from the electron recoils cause by background radioactivity. In contrast to 
other experiments, CRESST detectors allow the employment of a large variety of 
target materials, such as e.g. sapphire or tungsten. This allows a better sensitivity 
for detecting the WIMPs. Although they are already using sapphire, for the next 
project, CRESST II, it is planned to install a mass of about 10 kg, consisting 
of 33 modules of 300 g tungsten crystals. With 3 years of measurement time the 
experiment will be able to test a WIMP-nucleon cross section slightly smaller than 
the CDMS Soudan discussed above. 

The most sensitive de tecto r will be the GErmanium in liquid Nitrogen Under- 
ground Setup (GENIUS) DSS1 shown in Fig.^b (from Ref. 1170)1 . which will be able 
to test a WIMP-nucleon cross section as low as a » 10 -9 pb (see Fig. 0). Indeed 
such a sensitivity covers a large range of the parameter space of SUSY models 
with neutralinos as dark matter, as we will see in the next section. The GENIUS 
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DRIFT 




Fig. 10. Loft picture: DRIF detector sited on the Earth through the WIMP halo. Right picture: 
Schematic of TPC 

project is based on the idea to operate an array of 100 kg of Ge crystals directly 
in liquid nitrogen. The latter, which is very clean with respect to radiopurity, can 
act simultaneously as cooling medium and shield against external activities, using 
a sufficiently large tank ^12 metres in diameter at least. It has been shown using 
Monte Carlo simulations that with this approach the unwanted background is re- 
duced by three to four orders of magnitude. In order to demonstrate the feasibility 
of the GENIUS project a GENIUS Test-Facility (GENIUS-TF) shown in Fig. 
is under installation since March 2001 with 40 kg of detectors UlJJ It could test the 
DAMA evidence for dark matter by the annual modulation signature (similarly to 
the GEDEON project mentioned above) within about two years of measurement. 

Other experiments are e.g. the French project MACHe3 at the Joseph Fourier 
University, the PICASSO project at the University of Montreal, or the ORPHEUS 
project at the University of Bern, 

Efforts to build detectors sensitive to the directional dependence, i.e. recoil away 
from direction of Earth motion, has also been carried ou t. This i s an extension of 
the idea of annual modulation, but clearly more powerful 1 -^7^ 1 17^ 1 A detector sited 
on the Earth can look for the signal as shown in the left frame of Fig. (from 
Ref. 1170(1 . The detector will see the mean recoil direction rotate and back again 
over one day. Studies indicate that a WIMP signal could be identified with high 
confidence from as few as 100 detected WIMP-nucleus scattering events, compared 
with thousands for an annual modulation. 

Of course, to reconstruct such a three-dimensional direction is not simple, but a 
standard Time Projection Chamber (TPC) where the arrival time of the ionization 
signal is measured may be used for it. Although the energies of WIMP induced 
nuclear recoils are smaller than 100 keV, implying short ranges for the ionization 
tracks, this problem can be attacked using a gas such as Xenon or Argon at low 
pressure to extend these ranges to a few millimetres. The recoil energy and direc- 
tion can be measured by drifting the ionisation created along the recoil track to a 
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suitable readout plane, as shown in the right frame of Fig. 1101 (from Ref. I174( l. The 
low pressure and therefore the low target mass is not a problem, since as mentioned 
above only a small number of WIMP events are needed. Another advantage of this 
technique is that its resolution should allow virtually all Compton electrons to be 
rejected, since the electron range at a given total ionization is predicted to be many 
times longer than for a Xen on recoil. The Directional Recoil Identification From 
Tracks (DRIFT) projecfGZIlis precisely using this technique to detect dark matter. 
For the moment they are using a detector, DRIFT I, consisting of Carbon Disul- 
fide CS2 with 1 m 3 volume. It was installed at the Boulby mine in summer 2001. 
Although DRIFT-I is currently taking data, the next stage in the DRIFT project, 
DRIFT-II, is being planned. In order to achieve sensitivity to scattering cross sec- 
tions below 10 -6 pb, a scale up in mass is needed. DRIFT-II is thus proposed to 
have increased volume. 

Let us finally mention that links to web pages of the direct detection experiments 
mentioned in this section, and others experiments, can be found in the Dark Matter 
Portal http://lpsc.in2p3.fr/tep/fred/dm.html. The same web page can be used to 
look for indirect detection experiments. This is precisely the subject of our next 
discussion. 



4.2. Indirect detection 

There a re promising methods for the indirect detection of WIMPs in the halo 
| 175 | 17t> l p or example, WIMPs passing through the Sun and/or Earth may be 
slowed below escape velocity by elastic scattering. From subsequent scattering they 
will sink to the center and accumulate there, annihilating with another WIMP 
into quarks and leptons (see Fig. |2J), gauge bosons, etc. Decays of these annihi- 
lation products will produce energetic neutrinos that can be detected 11 . In par- 
ticular, in underground, underwater and under-ice experiments through the 'up- 
going' muons produced by their interactions in the rock, water, and ice respec- 
tively. Although some underground neutrino telescopes, such as e.g. Kamiokande 
or MACRO, have already reported null results, in the next years underwater (e.g. 
NESTOR, ANTARES) and under-ice (e.g. AMANDA) experiments with sizes of 
about 10 3 ' 4 m 2 (10 6 m 2 in the case ICECUBE, a proposed expansion of AMANDA) 
will have a much greater sensitivity. 

Another way of detecting WIMPs indirectly is through anomalous cosmic rays 
produced by their annihilations in the galactic halo 1 . Whereas positrons produced in 
this way are difficult to distinguish from the standard background of typical cosmic- 
ray positrons, there are other products of annihilation that can be distinguished. 
In particular, the most promising are cosmic-ray antiprotons and cosmic gamma 

h Another proposal for studying the WIMPs accumulated in the center of the Sun is through 
helioseismolog yllZZI 

'Detection of extragalactic cosmic gamma rays has also been analyzed recently l -^ 7? [ 
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rays. As discussed above, WIMPs annihilate into quarks, leptons, gauge bosons, etc. 
These will hadronize and produce, among other products, antiprotons. Unlike the 
usual cosmic-ray antiprotons, the flux of these does not decrease dramatically at 
energies less than a GeV. Space-based antimatter experiments, such as PAMELA 
and AMS, will study this possibility. On the other hand, two WIMPs can annihi- 
late into gammma rays through loop diagrams, and this might be detected using 
atmospheric Cherenkov telescopes (e.g. CELESTE, MAGIC) or space-based 7 ray 
detectors (e.g. AMS, GLAST). In particular through the observation of monochro- 
matic photons at an energy equal to the WIMP mass. No typical cosmic gamma 
rays are monochromatic at energies about 100 GeV. Another contribution of pho- 
tons is the continuum emission from all gamma rays produced by the cascade decays 
of other primarily annihilation products (e.g. from the 7r° decay originated by the 
fragmentation of quarks). Although the signal is less distinctive than the previous 
one, it has the advantage of being unsuppressed by loop effects, and therefore large 
flux of photons might be observed J . 

Finally, it is worth noticing that in some regions of SUGRA scenarios, exper- 
iments such as AMANDA, NESTOR, ANTARES, MAGIC, GLAST, A MS m ight 
be competitive with direct detection experiments such as CDMS Soudan LLlfll 



5. Theoretical models 

As discussed in the Introduction and Subsection 3.1.3, the leading candidate for 
WIMP is the neutralino, a particle predicted by SUSY extensions of the standard 
model. In this section we will review the known SUSY scenarios based on ncutralinos 
as dark matter candidates, in the context of the MSSM. In particular, we will 
discuss how big the cross section for the direct detection of these neutralinos can 
be. Obviously, this analysis is crucial in order to know the possibility of detecting 
the dark matter in current and projected experiments (see Fig.EJ. 

Since in the MSSM there are four neutralinos, Xi (i — L 2, 3, 4), the lightest, Xi> 
will be the dark matter candidate. The neutralinos are the physical superpositions 
of the fermionic partners of the neutral electroweak gauge bosons, called bino (B°) 
and wino (W3 ), and of the fermionic partners of the neutral Higgs bosons, called 
Higgsinos (H®, H®). The neutralino mass matrix with the conventions for gaugino 
and Higgsino masses in the Lagrangian, C = j^MjAjA,, + /j,H®H® + h.c, is 
given by 

/Mi -M z cos/3sin9 w M z sin p sin 9 W 

Mi Mz cos j3 cos 9-yy — Mz sin p cos 9w 

— M z cos P sin 9\y Mz cos P cos 9yy — /j, 

\ Mz sin P sin 9\y — Mz sin /3 cos 9y/ — \i 



J It is worth noticing that the EGRET telescope has identified a gamma-ray so urce at the galactic 
center that, apparently, has no simple explanation with standard proces ses \ ^J ^ [ Likewise, the 
HEAT balloon experiment has confirmed an excess of cosmic ray positrons | l?Q l 
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in the above basis (B° = -iX', W> 3 ° = -iX 3 , H°, Hf). Here tan/3 = (H°)/{H%) is 
the ratio of Higgs vacuum expectation values (VEVs), Mi and M 2 are the soft bino 
and wino masses respectively, and fj, is the Higgsino mass parameter determined by 
the minimization of the Higgs effective potential 

„ m 2 r, — rriir tan 2 B 1 „ 
K 2 = "WgLl " i"*' (13) 

Thus one parameterizes the gaugino and Higgsino content of the lightest neutralino 
according to 

Xl = N U S° + N 12 W° + N 13 H° d + N U H° . (14) 

It is commonly defined that Xi is mostly gaugino-like if P = |^Vn| 2 + |jVi 2 | 2 > 0.9, 
Higgsino-like if P < 0.1, and mixed otherwise. 

The relevant effective Lagrangian describing the elastic Y ?-nucleo n scattering in 
the MSSM has been examined exhaustively in the literature I I I It is given by 

£eff = a2iXl^l 5 X QiJ^Qi + auXXmi ■ (15) 

The Lagrangian is to be summed over the quark generations, and the subscript 
i refers to up-type quarks (i=l) and down-type quarks (i=2). The couplings 0:2,3 
can be found e.g. in Ref. 1281 with the sign conventions for Yukawa couplings in the 
Lagrangian, L = —X u H^uluji — XdH^d^dj^ — X e H^e.LeR + h.c. In particular, the 
coupling ct3i associated to the scalar (spin-independent) interaction include contri- 
butions from squark (g) exchange and CP-even light (h) and heav y (H ) neutral 
Higgs exchange, as illustrated in Fig. 0J and can be approximated as' 



23 



Xi Xi 

12 



99 m qi 
4mwB. 



I 1 \ — „ f£t . Cf 
ml 



SuNnDiQ —j - + 8 2i N n . , , 



H / \"*h ""H 



(16) 



where we are neglecting CP-violating phases k . Here, Su is Ni 3 (iVu), 82% is Nu 
(—N13), Bi is sin/3 (cos/3), C$ is sin a (cos a) with a is the Higgs mixing angle, 
and Di is cos a (—sin a) for up (down) type quarks, yi is the hypercharge defined 
by &i = T$i + yi/2. m qi are the masses of the quarks, and the masses mi and m 2 
correspond to the two squark mass eigenstates, and the ry's are the entries of the 
matrix diagonalizing the sfermion squared mass matrix which can be parameterized 
by an angle Of 

cos Of sin0/\ = ( Tin rj 12 \ 
- sinOf cosOf J \ 7721 7722/ 



k Analyses of the effects of the CP phases on the cross section can be found in Rcfs. 19 22 23 and 

121 
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The first two terms in Eq. I|16|) arise from interactions of the type g'qqB oc 
g' Nuqqxi, an d the last two from \ q H® d qq (note that gm q J (niwBi) = v^AgJ, 
g'H° d H° d B oc g'Ni S ,uNiiH° d XiXv Due to the relative size of the lightest Higgs 
mass to the squark masses, and the smallness of t]i2 to Tjn, the Higgs exchange term 
can dominate in 03 1 9 \ Note also that the heavy Higgs can make a significant 
contribution to the down- type-quark part of 013. The reason being that the latter 
is proportional to C d = cos 2 a whereas the light Higgs contribution is proportional 
to D\ — sin 2 a, and for a wide range of parameters one finds tana = 0(1/10). In 
addition, as we will discuss in Subsection 5.2.1, run may decrease significantly for 
large tan/?. 

The contribution of the scalar interaction to the x^-nucleus cross section, (J-^a_ Nl 
is given 

AM 2 

a m = ^[Zf p +(A-Z)f n } 2 , (17) 

where the reduced mass M r = M/vto^o/(M/v + m^o), with Mjv the mass of the 
nucleus, and the relevant contributions to f p are 

JjL - Ap)^1 1 2 f (p) ^3? n „s 

~ 2^ Jra + 27 Jtg2^ m > ^) 

q—u 7 d,s c,b,t 

where m p is the mass of the proton, m q the masses of the quarks, the parameters 
are defined by (p\m q qq\p) = m p f^ q \ while f$l = 1 - T, q =u,d,s frg - /» has a 
similar expression. Typical numerical values of the hadronic matrix elements f^ q n ^ 
can be found e.g. in Ref. 28J These are: 

f^l = 0.020 ± 0.004 , fjfy = 0.026 ± 0.005 , f$ = 0.118 ± 0.062 , 

= 0.014 ± 0.003 , ffij = 0.036 ± 0.008 , f$ = 0.118 ± 0.062 , (19) 

and we will use throughout this review the central values 1 . We see that f^J = fj?] 
and much larger than fx q for u and d quarks, and therefore f p and /„ are basically 
equal (note that the relevant part of the couplings a 3q are proportional to m q , so 
the fraction a 3q /m q does not become large for small m q ). Thus we can write™ 

a 3N ~ ^A 2 / 2 . (20) 

7T F 

It is worth noticing here that the spin-independent scattering adds coherently 
giving rise to a cross section proportional to the squared of the atomic weight, 
A. However, the axial-vector (spin-dependent) interaction, the one proportional to 
otn in Eq. I|15fl . which is nonzero only if the nucleus has a non- vanishing spin, is 

'Larger values of the hadronic matrix elements have also been considered in the literature. For 
analyses of the effect induced on the cross section by the present uncertainties in these values, see 
e.g. Refs, l2Hl and l29l A rough estimate can be obtained from Eq. 1231 below. 

m The contribution to this cross section from two-nucleon currents from pion exchange in the 
nucleus has recently been discussed in the last work of Ref. 11631 
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incoherent. Thus, for heavy targets, the scalar cross section is generically larger 
implying a%o_ N ~ 03./V. Recall that this is the case of the experiments discussed in 
Section 4. For example, Ge with a mass ~ 73 GeV and Nal with masses ~ 23 GeV 
for Na and ~ 127 GeV for I are used. In what follows we will concentrate on the 
scalar cross section 11 , and in particular on the one for protons 

rr ~rr 

CT 3p = °x\-p - ' ^ > 

where m r = m p m^o j (m p + m^a) « m p . In fact, this is the cross section shown in the 
experimental papers (see Fig. [JJ, and therefore the one that we will compute in the 
different theoretical models. Note that for a material with heavy nuclei, A rj 100, 
Mm w 100 GeV « m^o and therefore 

a x°-N ~ !0 8 a xl~v ■ ( 22 ) 

For c^o_ p m 10 -8 pb one recovers the rough estimate of Subsection 3.1.3, a^o_ N R* 
1 pb, for a particle with weak interactions. 

The predictions for the scalar neutralino-proton cross section, er^o_ p , are usually 
studied in the framework of SUGRA, and we will review them below. We will see 
that <J x "_ p w 10~ 8 pb can be obtained, but that smaller (or larger) values are 
also possible depending on the parameter space chosen (e.g. m, M, A, tan/3 in 
mSUGRA as we will discuss in Subsection 5.2.1), since very different values for m,h, 
run, Nn, N13, etc. can arise. Note in this sense that, from the above results, the 
neutralino-proton cross section can be approximated as 



a x° 1 -P 




99'fk 



rriw cos 



I 1 \ / sin 2 a 



NiiNi4s'macosa — j 5- + N11N13 



h m H/ v "'/, "'// 



2 

(23) 



For example, assuming the typical values N±i w 1, N13 w 0.1, tana w 1/10, 
tan/3 w 10, and m H w 100(500) GeV, one obtains cr^o_ p w 10~ 8 (10 -11 ) pb. 

Before analyzing in detail the cross section, it is worth remarking that in these 
analyses to reproduce the correct phenomenology is crucial. Thus, in the next Sub- 
section, the most recent experimental and astrophysical constraints which can affect 
this computation will be discussed. 



"Recent analyses of the spin-dependent part of the cross section can be found e.g. in Refs. l^3ll28l 
1341 1471 and 11841 In Ref. 118^1 it was pointed out that for targets with spin-non-zero nuclei, as e.g. 
73 Ge where spin = 9/2, it might be the spin-dependent interaction the one determining the lower 
bound for the direct detection rate, when the cross section of the scalar interaction drops below 
about 10~ 12 pb. 
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5.1. Experimental and astrophysical constraints 

We list here the most recent experimental and astrophysical results which are rel- 
evant when computing the neutralino-nucleon cross section. They give rise to im- 
portant constraints on the SUSY parameter space. 

(a) Higgs mass 

Whereas in the context of the standard model the negative direct search for the 
Higgs at the LEP2 collider implies a lower bound on its mass of about 114.1 GeV, 
the situation in SUSY scenarios is more involved. In particular, in the framework 
of mSUGRA (to be discussed in detail below), one obtains for the lightest CP-even 
Higgs nth > 114.1 GeV when tan/3 < 50, and nih > 91 GeV when tan/3 is larger 
HHSl Recall in this sense that <JsusY(e + e~ — > Zh) = sin 2 (a — 0) (7sA/(e + e~ — > Zh) 
with a the Higgs mixing angle in the neutral CP-even Higgs sector 1 °' ' Thus 
when the ZZh coupling, sin 2 (a — (3), which controls the detection of the lightest 
MSSM Higgs at LEP, is ~ 1 one recovers the standard model bound. However, 
when a suppression of sin 2 (a — (3) is obtained the bound will be smaller. Such a 
suppression occurs with tan/3 > 50. In this case, in some regions of the parameter 
space m 2 A = m 2 Hd + rr? H + 2/j? becomes small, rriA < 150 GeV, because the bottom 
Yukawa coupling entering in the renormalization group equation (RGE) for rn 2 Hd is 
large. And a small tua gives rise to sin 2 (a — /3) < 1, a s can be understood in terms 
of the relations between the angles a, /3 and tua USSI. 

In any case, let us remark that generically tan /3 is constrained to be tan /3 < 60, 
since otherwise several problems arise. For example, as mentioned above, the bottom 
Yukawa coupling is large and for tan/3 > 60, even for moderate values of m and 
M, mj ld becomes negative. As a consequence m 2 A becomes also negative unless a 
fine-tuning (in the sense that only certain combinations of m and M are possible) is 
carried out. In this review we study only cases with tan/3 < 50. Thus for mSUGRA 
we will always have sin 2 (a — (3) ~ 1. However, we will also be interested in relaxing 
the mSUGRA framework and therefore sin 2 (a — /3) must be computed in this case 
for all points of the parameter space, in order to know which bound for the lightest 
MSSM Higgs must be applied. For the latter one can use the plot sin 2 (a — /3) versus 
rrih shown in Ref. 11891 

Let us finally mention that a very convenient program to evaluate mt is the so- 
called FeynHiggsHSDJwhich contains the complete one-loop and dominant two-loop 
corrections. Higher-order corrections introduce an uncertainty of about 3 GeV in the 
result. In addition, there is a simplified version of the program, FeynHiggsFast. 
The value of nih obtained with this version is approximately 1 GeV below the 
one obtained using FeynHiggs. The figures shown throughout this review were 
obtained using FeynHiggsFast, and neglecting the uncertainty due to higher-order 
corrections. Recently, another program to study Higgs phenomenology has been 
constructed. See Ref. I191l for details. 

(b) Top mass 
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The central experimental value for the top mass, m t (pole) = 175 GeV, is used 
throughout this review. However, let us remark that a modification in this m ass by 
±1 GeV implies, basically, a modification also of ±1 GeV in the value of mh ' ®™ 

(c) Bottom and tau masses 

For the bottom mass the input mb{mb) = 4.25 GeV is used throughout this review, 
which, following the analysis of Ref. 11921 with a s (Mz) = 0.1185, corresponds to 
rrib(Mz) — 2.888 Ge V. In the evolution of the bottom mass the SUSY threshold 
corrections at MsuSY^^&re taken into account. These are known to be significant, 
specially for large values of tan/3. A similar analysis must be carried out for the 
tau mass, using as input m T {Mz) = 1.7463 GeV. 

(d) SUSY spectrum 

The present experimental lower bounds on SUSY masses coming from LEP and 
Tevatron must be imposed. In particular, using the low-energy relation from 
mSUGRA, Mi = |tan 2 #i^M 2 , one obtains for the lightest chargino mass the 
bound E31 m^± > 103 GeV. Likewise, one is also able to obtain the following 

bounds for sleptons masses ^123 TOg > 99 GeV, mp, > 96 GeV, m f > 87 GeV. 
Finally, for the masses of the sneutrino, the lightest stop, the rest of squarks, and 
gluinos, one can use the following bounds: > 50 GeV, > 95 GeV, rriq > 150 
GeV, rrig > 190 GeV. 

(e) b — > S7 

The measurements of B -> X sl decays at CLEO GUI and BELLE i ea d to 
bounds on the branching ratio b — > sj. In particular we impose throughout this 
review: 2 x 10~ 4 < BR(b — > s-f) < 4.1 x 10~ 4 . Let us mention tha t a ro utine to carry 
out this evaluation is provided e.g. by the program micrOMEGAsEIEl A description 
of this procedure can be found in Ref. 11991 Although the improvements of Ref. 200 
are not included in this routine, they are not so important for this review since, as 
discussed below, only p > will be considered. 

(!) g» - 2 

The new measurement of the anomalous magn etic moment of the muon, — (g^ — 
2)/2, in the E821 experiment at the BNlESH deviates by (33. 7± 11.2) x lO" 10 from 
the recent standard model calculation of Ref. 202 using e + e~ data. Assuming that 
the possible new physics is due to SUGRA, we will show in the figures throughout 
this review the constraint 11.3 x 10~ 10 < a^(SUGRA) < 56.1 x lO" 10 at the 2a 
level. This excludes the case /i < 0°. 

(g) B s -> (ti+ju- 



° However, it is worth noticing that this result for a M is in contradiction with the one obtained by 
using tau decay data (instead of e+e — ones) which only implies a deviation (9.4 ± 10.5) X 10 -10 
from the standard model calculation I222( In any case, concerning the value of fi, it is worth 
recalling that constraints coming from the b — » 57 process highly reduce the /i < parameter 
space. 
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The branc hing ratio for the B s — > /i + fi decay has been experimentally bounded by 
CDF 121121 with the result BR(B S ->■ fx + fi~) < 2 -6 x 10~ 6 . Although in the standard 
model this branching ratio is very small, of order 3 x 10 -9 , it might be in prin ciple 
significant in the SUSY case for large tan/3, due to Higgs (A) mediated decay ***™ . 
This issue has been analyzed recently in the context of dark matter. It seems that 
the current experimental constraint does not eliminate in fact any relevant part of 
the parameter space of the MSSM ESI. Analyses of other scenarios can be found in 
Refs.EUandEI 

(h) LSP 

As mentioned in the Introduction, the LSP, with mass of order GeV, is stable 
and therefore must be an electrically neutral (also with no strong interactions) 
particle, since otherwise it would bind to nuclei and would be detectable in the 
Earth as an exotic heavy isotope with abundanc e n/ n(r>roton) ~ 10 _10 (10 -6 ) in 
the case of strong (electromag netic ) interactions LUsHl This is not consistent with 
the experimental upper limits E0S1 n /n{proton) < 10~ 15 to 10~ 30 . Although the 
lightest neutralino, Xi> is the LSP in most of the parameter space of the MSSM, 
in some regions one of the staus, fx, can be lighter. Therefore, following the above 
discussion, these regions must be discarded. 

(i) Relic Xi density 

The preferred astrophysical bounds on the dark matter density (see Eq. (jSJ), 
0.1 < n, OM h 2 < 0.3, must be imposed on the theoretical computation of the relic 
Xi density, once fi^o = f2 DM is assumed 13 . For the sake of completeness, we also 
show in the figures below the bounds 0.094 < DM /i 2 < 0.129 deduced from the 
WMAP satellite (see Eq. As mentioned in footnote a, different cosmological 
scenarios give rise to different results in the computation of the relic density. We 
will consider throughout this review the standard mechanism of thermal production 
of neutralinos (see Section 3.1.3). 

Let us finally mention that there are programs I ^Q7|198 I to evaluate fl^o . A very 
convenient one used in the figures shown here is mi croMEGAsESEI. In this program 
the exact tre e-level cross sections for all possible annihilation^! and coan nihila tion 
| lUS I H)9 | l M channels are included in the code through a link to CompHEpGHEl anc [ 
accurate thermal average of them is used. Also, poles an d thr esholds are properly 
handled and one- loop QCD corrected Higgs decay widths E£M are used. The SUSY 
corrections included in the latest version of the code BED are no t implemented yet 



p If neutralinos do not constitute all the dark matter in the Universe, and other particles discussed 
in Section 3 also contribute, then Q~o < £lr>M- This may allow points in the parameter space which 
would be excluded by Q-o/i 2 < 0.1. However, notice that in this case the density of neutralinos in 
the galaxy would be p-o = 5 ' 0-3 GeV cm -3 with £ < 1, and therefore the results (event rates) 

would imply a larger value of the experimental cross sectionlMI Since the DAMA area is already 
difficult to reproduce in SUSY models, now the situation will be even worse. In Fig. |3 we should 
substitute o" WIMP _ nuc ; eori by £ ■ a w : M P _ n uc i eorh . 
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by micrOMEGAs. Fortunately, in our case, their effect is much smaller than that of 
the QCD corrections. Good agreement between micrOMEGAs and other independent 
computations of £l%a including Xi ~ t\ coannihilations can be found in Ref. 12101 

(j) UFB 

The constraints that arise from imposing the absence of charge and colour breaking 
minima can also be considered l^al As is well known, the presence of scalar fields 
with colour and electric charge in SUSY theories induces the possible existence 
of dangerous charge and colour breaking minima, which would make the standard 

jr'q 

vacuum unstable^. The presence of these instabilities may imply either that the 
corresponding model is inconsistent or that it requires non-trivial cosmology to jus- 
tify that the U niverse eventually fell in the phenomenologically realistic (but local) 
minimum There are two types of constraints: the ones arising from directions in 
the field-space along which the (tree-level) potential can become unbounded from 
below (UFB), and those arising from the existence of charge and color breaking 
(CCB) minima in the potential deeper than the standard minimum. By far, the 
most restrictive are the UFB bounds, and therefore these are the ones used in the 
figures shown here. There are three UFB directions, labelled as UFB-1, UFB-2, 
UFB-3 in Ref. 12121 It is worth mentioning here that in general the unboundedness 
is only true at tree-level since radiative corrections eventually raise the potential 
for large enough values of the fields, but still these minima can be deeper than 
the realistic one (i.e. the SUSY standard-model vacuum) and thus dangerous. The 
UFB-3 direction, which involves the scalar fields {H u , vr li , er-, e»-} with i ^ j and 
thus leads also to electric charge breaking, yields the strongest bound among all 
the UFB and CCB constraints. 

5.2. SUGRA predictions for the neutralino-nucleon cross section 

As can be deduced from the above discussion concerning the cross section in 
Eq. (|21(l . the values of the SUSY parameters, such as scalar masses or gaugino 
masses, are crucial in the analysis. Obviously, these values are associated with 
the mechanism of SUSY breaking. A pragmatic attitude to this issue is the ad- 
dition of explicit soft SUSY-breaking parameters of the appropriate size (of order 
10 2 — 10 3 GeV) in the Lagrangian and with appropriate flavour symmetries to 
avoid dangerous flavour-changing neutral currents (FCNC) transitions. The prob- 
lem with this pragmatic attitude is that, taken by itself, lacks any theoretical ex- 
planation. SUGRA theories provide an attractive context that can justify such 
a procedure. Indeed, if one considers the SUSY standard model and couples it 
to N = 1 SUGRA, the spontaneous breaking of local SUSY in a hidden sector 
generates explicit soft SUSY-breaking terms of the required form in the effective 
low-energy Lagrangian ESI If SUSY is broken at a scale As, the soft terms have a 
scale of order A s /Mpi anc k, and therefore one obtains the required size if SUSY is 
broken at an intermediate scale As ~ 10 11 GeV. 

Thus one usually considers the MSSM in the framework of SUGRA. Working in 
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this framework the soft parameters generated once SUSY is broken through gravi- 
tational interactions, i.e., gaugino masses, scalar masses, trilinear couplings and bi- 
linear couplings, are denoted at high energy by M a , m a , A a py, and B, respectively. 
Although in principle these are free parameters together with /i, when electroweak 
symmetry breaking is imposed /i 2 and B are determined. Finally, the renormal- 
ization group equations (RGEs) are used to derive low-energy SUSY parameters' 3 . 
Since the SUGRA framework still allows a large number of free parameters, in or- 
der to have predictive power, one usually assumes that the above soft parameters 
are universal at the GUT scale, Mqut ~ 2 x 10 16 GeV, providing also for an un- 
derstanding of FCNC suppression. As mentioned in the Introduction, this is the 
mSUGRA scenario, also called in the literature the Constrained MSSM (CMSSM). 
Although, apparently, this assumption about universality is arbitrary, it is worth 
noticing that interesting cl asses of SUGRA models give rise to this kind of universal 
soft SUSY-breaking terms^^l In addition, explicit string constructions with these 
universality properties can be found in some limits 1 * We will discuss them in 
some detail in Subsections 5.3 and 5.4. 

Let us recall that the full N=l Supergravity Lagrangian is specified in terms 
of three functions which depend on the scalar fields 4>m of the theory: the an- 
alytic gauge kinetic function f a {4>M), the real Kahler potential K(cj>M , 0m) an d 
the analytic superpotential W{4>m)- I n particular, f a determines the kinetic terms 
for the gauge fields and is related to the gauge couplings as Ref a — l/fl^, where 
the subindex a is associated with the different gauge groups of the theory, K de- 
termines the kinetic terms for the scalar fields, and W determines the Yukawa 
couplings. Thus, for example, the form of K that leads to canonical kinetic terms 
for the observable fields C a , namely K — ^ a C a C^, irres pectiv e of the SUSY- 
breaking mechanism gives rise to universal soft scalar masses ESI TOq = m, of the 
type imposed in mSUGRA. 

Below we will carry out the analysis of mSUGRA concerning dark matter. We 
will also discuss how the results are modified when the above assumptions are 
relaxed. In particular, we will allow an intermediate scale instead of the usual GUT 
one, and also non-universal soft scalar and gaugino masses. 

5.2.1. mSUGRA scenario with a GUT scale 

As can be deduced from the above discussion, in this scenario there are only four 
free parameters 1 : m, M, A, and tan/3. In addition, the sign of fi remains also 

q Let us remark that, given the convention used throughout this review for gaugino masses in 
the Lagrangian, C = i '}2 a M a \ a \ a + h.c, one has to use the (one-loop) RGEs obtained e.g. in 
Ref. l214l but with an opposite sign in the gaugino contributions to the RGE's of the A parameters. 
r In fact, scenarios with less free parameters are also possible. For instance, in the above example 
with canonical kinetic terms, if in addition we assume that Yukawa couplings and are constants, 
i.e. that they do not depend on hidden sector fields, one obtains B = A — m. Scenarios with this 
type of constraints between soft parameters have been analyzed in the context of dark matter in 
Ref.1551 
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Fig. 11. Running of the soft Higgs masses-squared with energy. 

undetermined (see Eq. (1131) ). Let us also remark that in this context the lightest 
neutralino is mainly bino* * '. To understand this result qualitatively, we show 
schematically in Fig. ^2 the well known evolution of m 2 ^ (and m 2 Hd neglecting 
bottom and tau Yukawa couplings) towards large and negative values with the scale. 
Since /i 2 given by Eq. i|13[) . for reasonable values of tan f3, can be approximated as, 

M 2 « ~m 2 Hu - \m% , (24) 

then it becomes also large (the effect of the one-loop corrections to the scalar 
potential can be minimized by evaluating fj, at the scale M$usy — y/ rn t 1 rn t 2 )- 
In particular, \/i\ becomes much larger than Mi and M^. Thus, as can be easily 
understood from Eqs. I|12|) and (|14|) . the lightest neutralino will be mainly gaugino, 
and in particular bino, since at low energy M\ = | tan 2 Oy^rM^ ~ 0.5A/2. We show 
this fact in the plot on the left frame of Fig. 1121 where for tan (3 = 10 the gaugino- 
Higgsino components-squared JV^ of the lightest neutralino as a function of its mass 
m^o are exhibited. Here we are using an example with m = 150 GeV and A = M. 
Note that M is essentially fixed for a given m^o. Clearly, Nn is extremely large 
and therefore P > 0.9. 

Now, using Eq. 121fl one can compute the cross section for different values of 
the parameters (for recent results see Refs. ITHin^ll2Tll2^1l2lll2lll2^ll2^ll2^1ISTl 

ESI Ei E£0 1301 111 111 113 ESI 113 mi mi EH EH E3 E3 EH1E11 E01 E3 El E3> 

As a consequence of the Xi being mainly bino, the predicted er^o_ p is well below 
the accessible experimental regions for low and moderate values of tan /3, since the 
scattering channels through Higgs exchange shown in Fig. El are not so important 
(recall that the Higgs-neutralino-neutralino couplings are proportional to -/V 13 and 
N14 as shown in Eq. ifTfljl ). In addition, the (tree-level) mass of the CP-odd Higgs, 
A, 

m A = m% t + m 2 Hu + 2p? , (25) 

will be large because /x 2 is large. Since the heaviest CP-even Higgs, H, is almost 
degenerate in mass with this, run will also be large producing a further suppression 
in the scattering channels. This fact is shown in Fig. 1131 where contours of c^o_ p 
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Fig. 12. Gaugino-Higgsino components-squared of the lightest neutralino as a function of its 
mass for the unification scale (left frame), Mi = 10 16 GeV, and for the intermediate scale (right 
frame), Mj = 10 11 GeV. 

in the parameter space (771, M) for tan/3 = 10, A = and /1 > s are plotted 
Note that we can deduce the value of the x? mass in the plots from the value of M, 
since m^o ~ Mi ~ 0.4 M. For the gluino mass we can also use the simple relation, 
nig ~ 2.5 M. 

As we can see in the figure, the experimental bounds discussed in Subsection 5.1 
are very important and exclude large regions of the parameter space. This is due to 
the combination of the Higgs mass bound with the g^ — 2 lower bound'. One obtains 
from the Higgs mass the lower bound M > 320 GeV, and from g^ — 1 the upper 
bound M < 440 GeV. These bounds imply for the neutralino, 128 < m^o < 176 
GeV, and for the gluino (and squarks) 800 < mg^ < 1100 GeV. The light shaded 
area in the figure shows the region allowed by the experimental bounds. There, the 
lower contour (double solid line) is obtained including also the constraint coming 
from the LSP bound, m^o < m^. For this area cr^o_ p w 10~ 9 pb. 

On the other hand, the astrophysical bounds 0.1 < fl^oh 2 < 0.3 must be im- 
posed in the computation. Let us recall that there are only four regions where the 
upper bound fl^oh 2 ~ l/cr ann < 0.3 can be satisfied. There is the bulk region at 

s Let us remark that the sign of the dominant contribution to the supersymmetric contribution 
to the — 2 is given by Mip-- As discussed in Subsection 5.1, we are taking this to be positive, 
and therefore we will only consider sign(M) = sign(fi). Now, due to the symmetry of the RGEs, 
the results for (— M, A, —fi) are identical to those for (M, — A,fi), and therefore, one can cover the 
whole (permitted) parameter space restricting to positive values for M and fi and allowing A to 
take positive and negative values. 

'Recall that we are using only the limit based on e+e~ analysis (see footnote n), otherwise the 
allowed region would extend towards the right hand side of the figure. 
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Fig. 13. Scalar neutralino-proton cross section cr-o_ J) in the parameter space of the mSUGRA 
scenario (m, M) for tan/3 = 10, A = and fi > 0. The dotted curves are contours of <r^o_ p . The 
region to the left of the near- vertical dashed line is excluded by the lower bound on the Higgs mass 
> 114.1 GeV. The region to the left of the near-vertical double dashed line is excluded by the 
lower bound on the chargino mass m_± > 103 GeV. The corner in the lower left shown also by a 

double dashed line is excluded by the LEP bound on the stau mass mf 1 > 87 GeV. The region 
bounded by dot-dashed lines is allowed by <? M — 2. The region to the left of the double dot-dashed 
line is excluded by b — » 57. From bottom to top, the solid lines are the upper bounds of the areas 
such as rrif < m-o (double solid), Q.~oh 2 < 0.1 and Q,~oh 2 < 0.3. The light shaded area is favored 

XI Xi Xi 

by all the phcnomcnological constraints, while the dark one fulfills in addition 0.1 < Q-oh 2 < 0.3 

Xi 

(the black region on top of this indicates the WMAP range 0.094 < fl- h 2 < 0.129). The ruled 

Xl 

region is excluded because of the charge and colour breaking constraint UFB-3. 



moderate M and m I I I There, since Xi is mainly bino, the annihilation chan- 
nels into leptons are important (see e.g. Fig. |2J, specially those from Ir exchange 
because they have the largest hypercharge. Moreover, the mass of Ir only receives 
very small contributions from gaugino loop diagrams. Here fl^oh 2 < 0.3 requires 
specific upper bounds on xi an d mas ses EIB Another region is the coannihi- 
lation region extending to larger M HUH. There the tail where the LSP is almost 
dege nerate with the NLSP, the stau, producing efficient coannihilations, is rescued 
LLiifll This is true even for large values of m^a . These two regions can clearly be seen 
in Fig. El There is also the focus-point region at m > 1 TeV not shown in this 
figure, and that will be discussed at the end of this subsection. Finally, there is the 
Higgs pole region extending to large M and m because rapid annihilation through a 
direct-channel pole HHH ^ 17 | 4U | wnen 2m^o ~ rriA,H (see e.g. Fig- El) - But this region 
occurs at large tan /3 and will be discussed below in Fig. 

Given the above bounds, we observe in Fig. ^| that the allowed (dark shaded) 
area becomes very small (extremely small if the recent WMAP data are taken into 
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Fig. 14. The same as in Fig. 1131 but for tan /3 

,2 



35 and 50. The white region at the bottom 



bounded by a solid line is excluded because m-^ becomes negative. 
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Scatter plot of the scalar neutralino-proton cross section <r^o_ p as a function of the 



Fig. 15 

neutralino mass m -o in the mSUGRA scenario, for tan f3 = 35 and 50, A = and fi > 0. The 
light grey dots correspond to points fulfilling all experimental constraints. The dark grey dots 
correspond to points fulfilling in addition 0.1 < fl-oh 2 < 0.3 (the black dots on top of these 

indicate those fulfilling the WMAP range 0.094 < n~ah 2 < 0.129). The circles indicate regions 
excluded by the UFB-3 constraint. The lines corresponding to the different experiments are as in 
Fig. E| 



account). Only the beginning of the tail mentioned above is rescued. In addition, the 
restrictions coming from the UFB-3 constraint exclude also this area. In conclusion, 
the results indicate that the whole parameter space for tan/3 = 10 is excluded on 
these grounds. This is also true for other values of A, as shown explicitly in Ref . 1631 
The n eutralino-pr oton cross section can be increased when the value of tan (3 is 
increased | 1& | 29 | 39|40 | Notice for instance that the contribution of the down-type 
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quark to the cross section is proportional to 1/ cos/3 (see e.g. Eq. (J22J)). In addi- 
tion, the bottom Yukawa coupling increases, and as a consequence Tn 2 Hd decreases, 
implying that m A , given by Eq. also decreases. Since, as mentioned above, 

ma ~ tua, this will also decrease significantly. Indeed, scattering channels through 
Higgs exchange are more important now and their contributions to the cross section 
will increase it. Thus, in principle, we can even enter in the DAMA region. However, 
the present experimental constraints exclude this possibility HZESI We show this 
fact for tan (3 — 35 and A = in the plot on the left frame of In principle, 

if we only impose the LEP lower bound m%± > 103 GeV, the cross section can be 
as large as &%a- p ~ 10~ 6 pb. However, at the end of the day, the other experimental 
bounds (Higgs mass, b — > 57, g^ — 2 upper bound) constrain the cross section to 
be c^i-p ~ 10~ 8 pb. The region allowed by the g^ — 2 lower bound is now larger 
and the cross section can be as low as Cj>°- P w 5 x 10~ 10 pb. These bounds im- 
ply 260 < M < 750 GeV, and therefore 104 < m^o < 300 GeV, 650 < m M < 1875 
GeV. Concerning the UFB-3 constraint, it is worth noticing that the larger tan /3 is, 
the larger the excluded region becomes. However, unlike the case tan/3 = 10, this is 
not sufficient to forbid the whole dark shaded area allowed also by the astrophysical 
bounds. In fact, one can check that with tan/3 > 20 one can always find regions 
where all constraints are fulfilled, depending on the value of A. For example, for 
tan (3 = 35 and A = M essentially the whole dark shaded area is allowed, whereas 
for A = — M, — 2M this is forbidden. A detailed analysis of this issue can be found 
in Ref . E3 

The above comments, concerning the cross section, can also be applied for very 
large values of tan /3, as e.g. tan f3 = 50 u . We show this in the plot on the right frame 
of Fig. O However, note that now, unlike the case tan/3 = 35, the whole dark 
shaded area allowed by experimental and astrophysical bounds is not constrained 
by the UFB-3. This is also true for other values of A. Let us finally recall that the 
region allowed by 0.1 < fl^oh 2 < 0.3 is larger because the CP-odd Higgs A becomes 
lighter as tan/3 increases, as discussed above. This allows the presence of resonances 
in the Higgs mediated annihilation channels, resulting in drastic reduction of the 
neutralino relic abundance. In the case of tan/3 = 50, the resonant effects in the 
annihilation channels are felt in the whole parameter space displayed in Fig. 1141 We 
can see as well, that the area of the parameter space where \i ~ ft coannihilations 
are relevant lead to values of Q^oh 2 < 0.1. 

We summarize the above results for tan/3 = 35,50, and A — 0, in Fig. 1151 
There, the values of cr%°-p allowed by all experimental constraints as a function 
of the neutralino mass m^o are shown. Dark grey dots correspond to those points 
having a relic neutralino density within the preferred range 0.1 < ilh 2 < 0.3. Given 
the narrow range of these points for the case tan/3 = 35, they overlap in the figure 

u For this very large values of tan/3 results become very sensitive to the precise values of mt and 
m,),. This has been discussed e.g. in Refs. 39 40 and 48 
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m. = 171 GeV , tan S = 10 , 11 > 

3500 ! , _ 1 




100 1000 1500 

mi/2 (GeV) 

Fig. 16. An expanded view of the M — m parameter plane showing the focus-point regions at large 
m for tan/3 = 10. In the shaded (mauve) region in the upper left corner, there are no solutions 
with proper electroweak symmetry breaking. Note that rat = 171 GeV has been chosen, in which 
case the focus-point region is at lower m than when mt = 175 GeV, as assumed in the other 
figures. The position of this region is very sensitive to mt ■ 



with those excluded by the UFB-3 constraint (shown with circles). We observe 
that, generically, the cross section and the neutralino mass are constrained to be 
(for any value of A) 5 x lO" 10 < a^o_ p < 3 x 10~ 8 pb and 120 < m^o < 320 GeV, 
respectively. 

Qualitatively similar results are obtained when dark matter is analyzed in the so- 
called focus-point supersymmetry scenario. Let us recall that this has been proposed 
as an alternative scenario in or der to avoid dangerous SUSY contributions to flavour 
and CP violating effects EEI The idea is to assume the existence of squarks and 
sleptons with masses which can be taken well above 1 TeV (of course this scenario 
rules out SUSY as an explanation of the possible deviation in the g M — 2 from the 
standard model prediction). It has also been argued that this situation produces 
no loss of naturalness. Notice that for m? » M 2 the electroweak scale given by 
Eq. J2U), M§/2 w —m 2 H ^ — /i 2 , can easily be obtained since rrijj becomes less 
negative. 

The implications of focus-point supersymmetry for neutralino dark matter have 
been considered in Ref. El In particular, it was pointed out that for m > 1 TeV 
the lightest neutralino is a gaugino-Higgsino mixture over much of parameter space. 
This can be understood from Eq. I|24|) since m 2 H becomes less negative for mo > 1 
TeV (recall Fig. Illfl . and therefore decreases. As a consequence, although as 
to increases the t-channel sfermion exchange process (see e.g. Fig. |2J) is more and 
more suppressed, the LSP gradually acquires a significant Higgsino component and 
other diagrams become unsuppressed. Thus the upper bound in Eq. @ is fulfilled. 
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This 'focus-point' region, which is adjacent to the boundary of the region where 
electroweak symmetry breaking is possible, is shown in Fig. Ilfil (from Ref. 49). 

Concerning the neutralino-proton cross section, scattering channels through 
Higgs exchange will increase it. However, one still needs large values for tan (3 (of or- 
der 50) in order to have very large cross sections, and then experimental constraints 
are very important. 

Obviously, in the mSUGRA scenario with a GUT scale that we have reviewed 
in this Subsection, where o~^o_ p ~ ^ x 10 8 P^i more sensitive detectors producing 
further data are needed. As discussed in Subsection 4.1.2, many dark matter de- 
tectors are being projected. Particularly interesting is the case of GENIUS, where 
values of the cross section as low as w 10~ 9 pb will be accessible, although this 
might not be sufficient depending on the values of the parameters (see Fig. I15|l . 

5.2.2. mSUGRA scenario with an intermediate scale 

The analysis of the cross section o~%o_ p carried out above in the context of 
mSUGRA, was performed assuming the unification scale Mqut ~ 10 16 GeV. How- 
ever, there are several interesting phenomenological arguments in favour of SUGRA 
scenarios with scales Mj w 10 10 ~ 14 G eV, such as to explain neutrino masses, the 
scale of axion physics, and others^^^. In addition, the string scale may be anywhere 
between the weak and the Planck scale, and explicit scenarios with intermediate 
scales may arise in the context of D-brane constructions from type I strings, as we 
will discuss in Subsection 5.3. Inspired by all these scenarios, to use the value of 
the initial scale Mj as a free parameter for the running of the soft terms is partic- 
ularly interesting. In fact, it was pointed out in Refs. 1331 1351 1411 that o-^o_ p is very 
sensitive to the variation of Mj for the running of the soft terms. For instance, by 
taking Mi = 10 10 ~ 12 GeV rather than M G ut, regions in the parameter space of 



Before trying to understand this result, let us discuss what we mean by an inter- 
mediate unification scale. Concerning this point two possible scenarios are schemat- 
ically shown in Fig. El for the example Mi = 10 11 GeV. In scenario (a) the gauge 
couplings are non universal, cti ^ a, and their values depend on the initial scale 
Mi chosen. As we will discuss in Subsection 5.3, a qualitatively similar scenario 
may arise in the context of type I string constructions if the gauge groups of the 
standard model come from different types of D-branes. Since different D-branes 
have associated different couplings, this implies the non universality of the gauge 
couplings. 

On the other hand, scenario (b) with gauge coupling unification at Mi, cti = a, 
can be obtained with the addition of extra fields in the massless spectrum. For the 
example of the figure these are doublets and singlets under the standard model 
gauge group. As discussed in Ref. ESJ the values of the gauge coupling constants 




xl-p 



is two orders of magnitude larger than for 
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Fig. 17. Running of the gauge couplings with energy, shown with solid lines, assuming (a) non 
universality and (b) universality of couplings at the initial scale Mj. For comparison the usual 
running of the MSSM couplings is also shown with dashed lines. 



at the intermediate scale are important in the computation of the cross section, 
and scenario (a) gives rise to larger cross sections than scenario (b). Let us then 
concentrate on the former. 

The fact that smaller scales imply a larger cr^o_ p can be explained with the 
variation in the value of /j, with Mr. One observes that, for tan /3 fixed, the smaller 
the initial scale for the running the smaller the numerator in the first piece of 
Eq. i|13|) becomes. This can be understood from the evolution of m? H with the 
scale (see Fig. Illll . Clearly, when the value of the initial scale is reduced the RGE 
running is shorter and, as a consequence, the negative contribution rr? H to pt 2 in 
Eq. H13JI becomes less important. Then, decreases and therefore the Higgsino 
composition of the lightest neutralino increases. Eventually, \fj,\ will be of the order 
of Mi , M2 and Xi WlU be a mixed Higgsino-gaugino state (see the plot on the frame 
of Fig. EJ- In addition, when |/z| decreases m A , given by Eq. jS5jl. also decreases. 
As mentioned in the previous Subsection when talking about increasing tan /?, H 
will decrease and therefore the scattering channels through Higgs exchange will 
increase the cross section. 

Let us also remark that, for the same value of the parameters, the lightest Higgs 
mass rrih decreases with respect to the GUT scale scenario. This is because the value 
of rrih depends on the value of the gluino mass M3. It increases when M3 increases 
at low energy. However, now the running is shorter and therefore M$ at low energy 
is smaller than in the GUT scenario. Although the latter may be welcome in order 
to obtain larger cross sections, it may also be dangerous when confronting with the 
experimental result concerning the Higgs mass. 

Concerning the value of the relic density, f^o is dramatically reduced with 
respect to the Mqut case. This is due to a combination of several factors: 1) 
The Higgsino-gaugino composition of \i allows a significant increase of the Xi 
annihilation cross section, due to channels with Higgs and gauge bosons in the 
final states; 2) The decrease of the mass of the pseudoscalar Higgs, along with the 
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value of the ^t-term, enables the presence of resonant annihilation channels even 
at tan/3 = 10; 3) T he masses of th e lightest chargino and stop are small enough 
to allow Xi~xt and Xi~ti coannihilations in some areas of the parameter 
space. Although the later is less relevant, we find some areas at tan/3 = 50 and 
A < 0. 

We show in Fig. [TBI the result for Mj = 10 11 GeV, with tan/3 = 10 and A = 
0, — 2M^l We choose A proportional to M becau se th is relation is particularly 
interesting, arising naturally in several string models However our conclusions 
will be independent on this assumption. For example, if we choose to do the plots for 
different constant values of A, a very common procedure in pure SUGRA analyses, 
the results will be qualitatively similar. The case A = can be compared with the 
one in Fig. El where Mqut is used. Now the relation m^o ~ 0.4 M does not hold, 
and one has m^o > 0AM. In any case, m^o < Mi since the bino-Higgsino mixing 
is significant in this case. Clearly, for the same values of the parameters, larger 
cross sections can be obtained with the intermediate scale. It is worth noticing that 
even with this moderate value of tan/3, tan/3 = 10, there are regions where the 
cross section enters in the DAMA area, cr^o_ p w 10 -6 pb. However, for A = 
the whole parameter space is forbidden due to the combination of the Higgs mass 
bound with the g M — 2 lower bound (for A = M the situation is similar). We have 
checked explicitly that 114.1 GeV is the correct lower bound to be used concerning 
the Higgs mass, since generically sin 2 (a — (3) ~ 1 for the intermediate scale. Notice 
also that now, for A = (and also for A = M), Q^oti 2 is smaller than 0.1 in most 
of the parameter space. Only tiny regions bounded by solid lines in the figure, and 
therefore with 0.1 < f^o/i 2 < 0.3, can be found. 

On the other hand, for A = —2M (and also for A = —M) there are small 



February 1, 2008 22:9 WSPC/INSTRUCTION FILE 

darkijmpa21+refs+refs+refs 



Dark Matter Detection in the Light of Recent Experimental Results 45 
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Fig. 19. The same as in Fig. 1151 but for the intermediate scale Mj = 10 11 GeV, with tan/3 = 
10, 35 and A = -2M. 



regions where the irih and g M — 2 bounds are compatible, but finally the constraint 
nih > 114.1 GeV implies that the allowed cross sections do not enter in the DAMA 
area. Although now larger regions with 0.1 < fl^oh 2 < 0.3 are present, for A = —M 
these are not compatible with the experimental bounds. 

One also finds that the regions excluded by the UFB-3 constraint are much 
smaller than in those cases where the initial scale is the GUT one. For example for 
A = (and A = M) no region is excluded (see however Fig. ^] for the GUT case). 
For A = —M the region excluded by the UFB-3 is smaller than the one forbidden 
by the LSP bound. We have to go to A = —2M to have it larger. 

In the plot on the left frame of Fi g. 1191 we summarize the above results for 
tan /3 = 10, concerning the cross section'""', showing the values of cr^o_ p allowed by 
all experimental constraints as a function of the ncutralino mass m^o , for A — —2M. 
Only in this case there are dark grey dots corresponding to points having a relic 
neutralino density within the preferred range 0.1 < Q^oh 2 < 0.3. Given the narrow 
range of these points, they overlap in the figure with those excluded by the UFB-3 
constraint. They correspond to 150 < m$o < 250 GeV, and e.g. the mass of the 
lightest stop, t-y, is between 250 and 340 GeV. 

Qualitatively, similar results are obtained for larger values of tan [3. For example, 
for tan j3 = 35 only for A = —2M we obtain points allowed by all experimental and 
astrophysical constraints, and this is shown in the plot on the right frame of Fig. 1191 
The upper bound in the cross section is because of the b — > s"f process. Points within 
the preferred astrophysical range correspond to cr^o_ p < 10~ 8 pb. Now, there are 
two allowed regions with 275 < < 325 GeV and 370 < m.o < 420 GeV. For 
these, 570 < m t - < 720 GeV. Let us finally mention that in the case of tan /3 — 
50, for A = —M there are points allowed by all experimental and astrophysical 
constraints, and cr^o_ p < 10~ 7 pb. 

It is worth noticing that in this analysis gaugino mass universality has been 
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assumed at the high-energy scale, although in this scenario gauge couplings do not 
unify. This situation is in principle possible in generic SUSY models, however it 
is not so natural in SUSY models from SUGRA where gaugino masses and gauge 
couplings are related through the gauge kinetic function. Since an explicit string 
construction with nonuniversal gauge couplings and gaugino masses will be analyzed 
in detail in Subsection 5.3, we have chosen to simplify the discussion here assuming 
gaugino mass universality. 

Summarizing, when an intermediate scale is considered in mSUGRA, although 
the cross section increases significantly the experimental bounds impose o~^o_ p < 4x 
10~ 7 pb. And, in fact, at the end of the day, the preferred astrophysical range for 
the relic neutralino density, 0.1 < fi^o/i 2 < 0.3, imposes cr^o_ p < 10~ 7 pb. Clearly, 
present experiments are not still sufficient, and more sensitive detectors producing 
further data are needed, as in the case of a GUT scale. 



5.2.3. SUGRA scenario with non-universal soft terms 

The gene ral si tuation for the soft parameters in SUGRA is to have a non-universal 
structure EEI Yor the case of the observable scalar masses this is due to the non- 
universal couplings in the Kahler potential between the hidden sector fields breaking 
SUSY and the observable sector fields. For example, K = ^ Q K a (h m , h^) C a C^, 
with K a a function of the hidden-sector fields h m , will produce non-universal scalar 
masses m a ^ mp if K a ^ Kp. For the case of the gaugino masses this is due to the 
non-universality of the gauge kinetic functions associated with the different gauge 
groups f a (h m ). For example, non-universal gaugino masses are obtained if the f a 
have a different dependence on the hidden sector fields. We will see in Subsection 5.3 
that gener al st ring constructions whose low-energy limit is SUGRA, exhibit these 
properties I^EI 

It was shown in the literature that the non-universality of the soft parameters 
allows to increase the neutralino-proton cross section with respect to the universal 
case. This can be carried out with non-universal scalar masses and/or gaugino 
masses We will concentrate on this possibility here. 



(i) Non-universal scalar masses 



Let us analyse a SUGRA scenario with GUT scale and non-universal soft scalar 
masses. In fact, non-universality in the Higgs sector, concerning dark matter, was 
first studied in Refs. 02171 12211 P*8 Subsequently, non-universality in the sfermion 
sector was added in the analysisE232EHI Analyses of the dark matter cross section 
using generic non-universal soft masses were carried out in Refs. 1181 1211 1291 1341 1451 
EOIESIESI whereas in Refs. ESIESI EH ED BHHH SU(5) or 50(10) GUT relations 
were used. In the light of the recent experimental results, one important consequence 
of the non-universality is that the cross section can be increased in some regions of 
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the parameter space. 

Let us then parameterized this non-universality in the Higgs sector as follows: 

m 2 Hd =m 2 (l + 5i) , m£ n = m 2 (l + 5 2 ) . (26) 

Concerning squarks and sleptons, in order to avoid potential problems with FCNC, 
one can assume that the first two generations have a common scalar mass m at 
Mgut, and that non-universalities are allowed only for the third generation: 

m QL = r " 2 ( 1 +<k) > m « H = m2 ( 1 + ^4), 
■m 2 eR = m 2 (l + 6 5 ) , mi B = m 2 (l +S 6 ) , 
m 2 LL = m 2 (l + <5 7 ) , (27) 

where Ql — Ll — Ur = £r and en = tr. Note that whereas 

8i > — 1, i = 3,..., 7, in order to avoid an UFB direction breaking charge and 
colour, <5i 2 < — 1 is possible as long as the conditions m 2 — m 2 H + fi 2 > 0, 
m 2 = m lr + M 2 > are fulfilled. 

As discussed for intermediate scales in Subsection 5.2.2, an important factor 
in order to increase the cross section consists in reducing the value of This 
value is determined by condition (|24|l and can be significantly reduced for some 
choices of the S's. We can have a qualitative understanding of the effects of the 
S's on fx from the following. First, when m 2 H at Mqut increases its negative low- 
energy contribution to Eq. (|24|l becomes less important. Second, when rriQ L and 
m 2 R at Mqut decrease, due to their contribution proportional to the top Yukawa 
coupling in the RGE of rr? H , the negative contribution of the latter to /i 2 is again 
less important. Thus one can deduce that /i 2 will be reduced (and hence c^o_ p 
increased) by choosing £3^ < and 62 > 0. In fact non-universalities in the Higgs 
sector give the most important effect, and including the one in the sfermion sector 
the cross section only increases slightly. Thus in what follows we will take Si = 0, 
i = 3, ...,7. 

Concerning the value of the relic density, fi^o is affected due to the increase of 
the Higgsino components of \i with respect to the dominant bino component of the 
universal case. The change in fi also determines the presence of the Higgs mediated 
resonant channels. In contrast to Subsection 5.2.2, the most relevant coannihila- 
tion scenarios are Xi~Ti, m particular Xi~Xi coannihilations are only sizeable for 
tan/3 = 35 in Fig. |201 (see the discussion below), when the /Lt-parameter becomes 
small. However, even in this case the area inside the WMAP bounds corresponds 
to neutralino annihilations, which are enhanced due to enlargement of its Higgsino 
components. 

On the other hand, there is another relevant way of increasing the cross section 
using the non-universalities of the Higgs sector. Note that decreasing rrijj , i.e. 
choosing Si < 0, leads to a decrease in m 2 A given by Eq. i|25|). and therefore in 
the mass of the heaviest Higgs H (on the contrary, the lightest Higgs mass, nth, is 
almost unaltered, it only decreases less than 1%). This produces an increase in the 
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Fig. 20. The same as in Fig. [Q] but for the non-universal case a) <5i = 0, &2 = 1, discussed in 
Eq. with tan/3 = 35, 50 and A = 0. 
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Fig. 21. The same as in Fig. [El but for the non-universal case a) <5i = 0, 82 = 1, discussed in 
Eq. with tan/3 = 35, 50 and A = 0. 



cross section v . 

Thus we will see that, unlike the universal scenario in Subsection 5.2.1., with 
non-universalities is possible to obtain large values of the cross section, and even 
some points enter in the DAMA area fulfilling all constraints. Let us analyze three 



v This effect might also be important when non-universal gaugino masses are taken into account. 
The contribution of M3 proportional to the bottom Yukawa coupling in the RGE of m \ d wm d° 

this smaller if M3 is large 
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Fig. 22. The same as in Fig. 1151 but for the non- universal case b) <5i = — 1, 82 = 0, discussed ir 
Eq. with tan/3 = 35, 50 and A = 0. 
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Fig. 23. The same as in Figs. El but for the non-universal case c) 8\ = —1, 82 = 1, discussed in 
Eq. with tan/3 = 35, 50 and A = 0. 



representative cases <^ with 



a) Si = 


. s 2 


= 1 


b) h = -1 


, s 2 


= 


c) Sx = -1 


, s 2 


= 1 



Clearly, the above discussion about decreasing /i 2 applies well to case a), where 
the variation in m 2 H through 5 2 is relevant. This is shown in Fig. I2UI for tan/3 = 
35, 50 and A = 0, which can be compared with Fig.^J Note that now, for tan/3 = 
35, there is an important area in the upper left where /i 2 becomes negative due to 
the increasing in S 2 with respect to the universal case. A larger area is forbidden for 
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large values of tan (3, as e.g. tan j3 = 50, but now because m 2 A becomes negative. This 
is similar to what occurs in the universal scenario when tan (3 > 60, as discussed 
in Subsection 5.1. Notice that from Eq. (|24|l one can write m 2 A in Eq. (|25J) as 
m A w rrijj — rrijj^ — M|. Since rr? Hu at Mqut increases its negative low-energy 
contribution becomes less important. In addition, the bottom Yukawa coupling is 
large and the m 2 H becomes negative. As a consequence m A < 0. 

For tan/3 = 35, although the cross section increases with respect to the uni- 
versal case, and is generically above the GENIUS lower limit, the present exper- 
imental constraints exclude points entering in the DAMA area. This can be seen 
more clearly comparing Figs. I21landll5l Notice also that the astrophysical bounds 
0.1 < il^oh 2 < 0.3 imply cr<>o_ p w 10~ 8 pb. On the contrary, for tan/3 = 50 there 
are points entering in the DAMA area, and even part of them fulfil the astrophys- 
ical bounds. We have checked that for A = M the figures are similar, although 
no points enter in the DAMA area, even for tan/3 = 50. On the other hand, the 
region forbidden by the LSP bound is larger than the one forbidden by the UFB-3 
constraint. 

We have also checked that larger values of 82, as e.g. 82 — 1.5, give rise to 
similar figures. For small values, 82 <^ 0.2 is sufficient to enter in DAMA fulfilling 
the experimental bounds with tan/3 = 50. In fact, e.g., for 62 — 0.5,0.75 one also 
gets many points entering in DAMA as for 82 = 1, however, they do not fulfil the 
astrophysical bounds. For the latter one needs 82 > 0.85. 

Let us finally remark that although sin 2 (a — (3) is close to 1 in most of the 
points, some of them can have smaller values when tan (3 = 50. As discussed in 
Subsection 5.1, these must be points with small values for tua, and in fact in this 
case are those close to the region with m 2 A < 0. The same situation occurs for the 
other cases studied below. Thus, according to our discussion in Su bsection 5.1, we 
use for these points the appropriate bound on the Higgs mass^^. In particular, in 
Fig. |21 the light grey dots above the CDMS line correspond to these points. 

For case b) the cross section increases also substantially with respect to the 
universal case. Now 82 is taken vanishing and therefore the value of fj, is essentially 
not modified with respect to the universal case. However, taking 8± = — 1 produces 
an increase in the cross section through the decrease in m 2 A , as discussed previously. 
As shown explicitly in Fig. |22 for tan (3 = 35 and ^4 = 0, there are points in the 
DAMA region. All of them correspond to sin 2 (a — 0) not close to 1, and therefore 
with an experimental bound on the Higgs mass smaller than 114.1 GeV. All points 
with sin 2 (o! — (3) ~ 1 have cr^o_ p < 6 x 10~ 7 pb. Note that points with large values 
of the cross section fulfil in this case the astrophysical bound fi-^o/i 2 > 0.1. Large 
values of m reduce the resonant effects produced by the Higgs A, and are sufficient 
to place the relic abundance inside the bounds. For tan/3 = 50, similarly to case a), 
there are points entering in the DAMA area, and part of them fulfil the astrophysical 
bounds. Those above the ZEPLIN line correspond to sin 2 (a — (3) not close to 1. For 
A = M the figures as similar. 
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We have checked that smaller values of 5\, as e.g. Si = —1.5, —2, give also rise to 
similar figures. For larger values, Si < — 0.4 is sufficient to enter in DAM A fulfilling 
the experimental and astrophysical bounds with tan/3 = 50. 

Finally, given the above situation concerning the enhancement of the neutralino- 
proton cross section for a) and b), it is clear that the combination of both cases 
might be interesting. This is carried out in case c) where we take S± = — 1 and 
#2 = 1. As shown in Fig. 1231 cross sections as large as o^o-p > 10~ 6 pb, entering in 
DAMA and fulfilling all experimental and astrophysical bounds, can be obtained for 
tan (3 = 35, 50 and A — 0. Those above the ZEPLIN line correspond to sin 2 (a — /3) 
not close to 1. On the other hand, for A — M and tan/3 = 35 no points with the 
correct relic density enter in DAMA. For other cases the results are similar. For 
example, if we consider Si = —0.5 and S2 = 1, one obtains for tan/3 = 35 points 
entering in DAMA but with tt^oh 2 < 0.01. For tan (3 = 50 a similar plot to the one 
in Fig. ESI is obtained. 

Concerning the restrictions coming from the UFB-3 constraint, we can see in 
Fig. [20| that these are slightly less important than in the universal scenario (see 
Fig. I14fl . Of course, this is not a general result, and different choices of the 5's can 
modify the situation. For example, for the same case as in Fig. [201 with tan/3 = 35, 
but using the opposite choice for the sign of the S parameters, not only the cross 
section is smaller, cr^o_ p < 10~ 8 pb, but also the UFB-3 constraint is very restric- 
tive, forbidding all points which are allowed by the experimental and astrophysical 
constraints. 

In summary, when non-universal scalars are allowed in SUGRA, for some special 
choices of the non-universality, the cross section can be increased a lot with respect 
to the universal scenario. It is even possible, for some particular values of the 
parameters, to find points allowed by all experimental and astrophysical constraints 
with cr^o_ p 10~ 6 pb, and therefore inside the DAMA area. Note however that 
these points would be basically excluded by the other underground experiments. In 
any case, the interesting result is that large regions accessible for future experiments 
are present. 

(ii) Non-universal gaugino masses 

The effects of the non-universality of gaugino masses on the dark matter relic 
density in SUGRA scenarios were studied in detail in Refs. \TM ITU71 l2l?51 ITTT1 
Analyses of the neutralino-proton cross section using SU(5) GUT relations for 
gaugino masses were carried out in Refs. |3U1 E21 whereas in Refs. LH3 EH1 
1561 E3I generic non-universal soft masses were used. It was also realized that the 
non-universality in the gaugino masses can increase the cross section. 

Let us parameterize this non-universality at Mqjjt as follows: 



Mi = M(l + S[) , M 2 = M (1 + 5' 2 ) , M 3 = M(l + 5' 3 ) , 



(29) 
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Fig. 24. The same as in Fig. 1151 but for the case discussed in Eq. 1291 with non- universal soft 
gaugino masses, 8[ 2 = 0, <5g = —0.5. 

where Mi^,3 are the bino, wino and gluino masses, respectively. Let us discuss now 
which values of the parameters are interesting in order to increase the cross section 
with respect to the universal case 8^ = 0. In this sense, it is worth noticing that 
M3 appears in the RGEs of squark masses, so e.g. their contribution proportional 
to the top Yukawa coupling in the RGE of vr? H will do this less negative if M3 is 
small. Although this effect increases the cross section, it is also worth noticing that 
small values of M3 also lead to an important decrease in the Higgs mass (we have 
checked that generically the value of sin 2 (a — f3) is very close to 1, and therefore 
we are using the lower bound rrih — 114.1 GeV as in the mSUGRA scenario). In 
addition, b — > sj and — 2 constraints are also relevant. 

Summarizing, although the cross section increases with respect to the universal 
case, the present experimental c onst raints exclude points entering in the DAMA 
region. This is shown in Fig. ^ for tan/? = 35,50 and A = 0, using S[ 2 = 
0, 5' 3 — —0.5, where one can see that there are points allowed by all experimental 
and astrophysical constraints, but they correspond to cr^.o_ p < 10 -7 pb. 

Finally, as in the previous case with non-universal scalars, increasing the cross 
section through values at low energy of rn 2 H less negatives implies less important 
UFB constraints. Now these are not very relevant, and in fact they correspond to 
the UFB-1 ones. 

5.2.4. Effective MSSM scenario 

In Refs. 1161 1251 1361 1371 1551 the authors considered that the uncertainties involved 
in SUGRA scenarios (such as e.g. the values of the soft parameters, i.e. universal- 
ity versus non-universality, or the choice of the scale, i.e. Mqjjt versus Mf) were 
problematic enough as to work better with a phenomenological SUSY model whose 
parameters are defined directly at the electroweak scale. This effective scheme of 



February 1, 2008 22:9 WSPC/INSTRUCTION FILE 

darkijmpa21+refs+refs+refs 



Dark Matter Detection in the Light of Recent Experimental Results 53 

the MSSM was denoted by effMSSM in Ref. 03 For example, in this work it was 
imposed for simplicity a set of assumptions at the electroweak scale: a) all trilin- 
ear parameters are set to zero except those of the third family, which are unified 
to a common value A; b) all squark soft-mass parameters are taken degenerate: 
niq i = nig-, c) all slepton soft-mass parameters are taken degenerate: mj = mf, d) 
the gaugino masses, Mi and are assumed to be linked by the GUT relation 
Mi = | tan 2 6\yM2 (this assumption was relaxed in the fifth and fourth works of 
Ref. 1871 obtaining relic neutralinos significantly lighter than those commonly con- 
sidered, m^o > 50 GeV). As a consequence, the SUSY parameter space consists of 
seven independent parameters, chosen them to be: M2, /x, tan/3, uia, rriq, m~ v A. 

Obviously, the scenarios studied in the previous Subsections, where some bound- 
ary conditions at the high scale are imposed and then the RGEs are used in order to 
obtain the low-energy spectrum, should be a subset of the general effMSSM scenario 
(i.e. without any assumption). In the effMSSM with the above assumptions, after 
imposing the experimental and astrophysical constraints, points with cross sections 
entering in the DAMA area can be obtained. This is similar to the result obtained 
when non-universal soft scalar masses were considered in the previous Subsection. 
Note in this sense that varying /x and ttia arbitrarily corresponds to vary mu u and 

m H d - 



5.3. Superstring predictions for the neutralino-nucleon cross 
section 

Although the standard model provides a correct description of the observable world, 
there exist, however, strong indications that it is just an effective theory at low 
energy of some fundamental one. The only candidate for such a theory is, nowadays, 
the string theory, which have the potential to unify the strong and electroweak 
interactions with gravitation in a consistent way. 

In the late eighties, working in the context of the (perturbative) Eg x Eg het- 
erotic superstring, a number of interesting four-dimensional vacua with pa rticle 
content not far from that of the supersymmetric standard model were found E2E1 
Until recently, it was thought that this was the only way in order to construct 
realistic superstring models. However, in the late nineties, it was discovered that 
explicit models with realistic prope rties can also be constructed using D-brane con- 
figurations from type I string vacua^23 We will review below these two superstring 
constructions concerning dark matter detection. 

Such constructions have a natural hidden sector built-in: the complex dilaton 
field S and the complex moduli fields T,;. These gauge singlet fields are generically 
present in four-dimensional string models: the dilaton arises from the gravitational 
sector of the theory and the moduli parameterize the size and shape of the com- 
pactified variety. Therefore the auxiliary fields of those multiplets can be the seed 
of SUSY breaking, solving the arbitrariness of SUGRA where the hidden sector is 
not constrained. In addition, in superstrings, K and / can be computed explicitly 



February 1, 2008 22:9 WSPC/INSTRUCTION FILE 

darkijmpa21+refs+refs+refs 



54 Carlos Munoz 

leading to interesting predictions for the soft parameters, and therefore for the value 
of the neutralino-proton cross section. 

5.3.1. Eg x Eg heterotic superstring constructions 

For any four-dimensional construction coming from the perturbative 10-dimensional 
heterotic superstring, the tree-level gauge kinetic function is independent of the 
moduli sector and is simply given by 

fa = k a S , (30) 

where usually one takes k 3 = k 2 = | k\ = 1. In any case, the values k a are irrelevant 
for the tree-level computation since they do not contribute to the soft parameters. 

On the other hand, the Kahler potential has been computed for several compact- 
ification schemes. This is for example the case of 6-dimensional Abelian orbifolds, 
where three moduli Tj are generically present. For this class of models the Kahler 
potential has the form 

K = - log(S + S*) - l°g( T * + T?) + ]T IC^UiiTi + T*)< . (31) 

i a 

Here n l a are (zero or negative) fractional numbers usually called "modular weights" 
of the matter fields C a . 

It is important to know what fields, either S or Tj, play the predominant role in 
the process of SUSY breaking. This will have relevant consequences in determining 
the pattern of soft parameters, and therefore the spectrum of physical particles. 
That is why it is very useful to introduce the f ollowing; parameterization for the 
VEVs of dilaton and moduli auxiliary fields EEI 

F s = V3{S + S*)m 3/2 sin 9 , 

F* = V3(T, + T*)m m cos 9 4 , (32) 

where i = 1, 2, 3 labels the three complex compact dimensions, m 3 / 2 is the gravitino 
mass, and the angle 8 and the 8.; with J2i l@j| 2 = 1) j us * parameterize the direction 
of the goldstino in the S, Ti field space. 

Using this parameterizatio n an d eqs. H30(l and (|31(l one obtains the following 
results for the soft parameters EEfl 



M a = V3; 



sm ( 



ml =m 2 — 2 ° V " ,q2 



3/2 ^l + 3cos 2 ^<e 2 ^ 

A a01 = -V3m 3/2 ^sin 9 + cos © 4 [l + < +n} 3 + nl / - (T 4 + T*)d % log Y a/3 ^ . 

(33) 

Although in the case of the A-parameter an explicit TVdependence may appear 
in the term proportional to di logY^g 7 , where Y a /3j(Ti) are the Yukawa couplings, 
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it disappears in several interesting cases, and we will only consider this possibility 
here. Using the above information, one can analyze the structure of soft parameters 
available in Abelian orbifolds. 

In the dilaton-dominated case (sin# = 1) the soft parameters are universal, and 
fulfil 

m = m 3/2 , M = V3m 3/2 , A = —M . (34) 

Of course, they are a subset of the parameter space of mSUGRA, and therefore one 
should expect small dark matter cross sections, as discussed in Subsection 5.2.1. 
However, in general, the soft parameters show a lack of universality due to the 
modular weight dependence (see scalar masses and trilinear parameters in Eq. (|33|l l. 
For example, assuming an overall modulus, i.e. T = Tj and Oi = 1/v3j one obtains 

m a — m 3/2 (l + n a C °S 2 <?) > 

Aa/3^ = — v3m3/2 sin# — m.3/2 cos# (3 + n a + np + n 7 ) , (35) 

where we have defined the overall modular weights n a — n l a , and in the case 
of Z n Abelian orbifolds they can take the values -1,-2,-3,-4,-5. Fields belonging to 
the untwisted sector of the orbifold have n a = —1. Fields in the twisted sector but 
without oscillators have usually modular weight -2 and those with oscillators have 
n a < —3. Note that e.g. for scalars in the untwisted sector 

m l + m } + m 7 = Af2 - A = ~ M ■ ( 36 ) 

These type of relations between fields in the same Yukawa coupling, for instance 
m QL + m u R + TO ff„ = were applied in Ref. 1^51 to the analysis of dark matter 
detection, with the final result that the neutralino-proton cross section is very 
similar to the one of the mSUGRA scenario. 

On the other hand, the apparent success of the joining of gauge coupling con- 
stants at Mqut ~ 2 x 10 16 GeV in the MSSM is not automatic in the heterotic su- 
perstring, where the natural unification scale is Mh — \p$Mpi anc k, with aw 1/24 
the gauge coupling. Thus unification takes place at energies around a factor 12 
smaller than expected in the heterotic superstring. This problem might be solved 
with the presence of large string threshold corrections which explain the mismatch 
between Mqut and Mh- In a sense, what would happen is that the gauge coupling 
constants will cross at Mqut and diverge towards different values at Mh- These 
different values appear due to large one-loop stringy threshold corrections. It was 
found in Ref. 2.28 that these corrections can be obtained for restricted values of the 
modular weights of the fields. In fact, assuming flavour independence, one finds that 
the simplest possibility corresponds to taking the following values for the standard 
model fields: 

tiq l = n dR = -1 , n UR = -2 , ul l = n eR = -3 , n Hu + n Hd = -5, -4 , (37) 
where e.g. ur denotes the three family squarks ur, cr, tn- 
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The above scenario was studied in Ref. 229 for uh u = — 3,n# d = —2, paying 
special attention to the calculation of the relic neutralino density. Subsequently, in 
Ref. El the authors estimate the direct detection rate for the same scenario. The 
conclusion was that this is small and large scale detectors are needed to discover 
the neutralino. 

Although, apparently, gaugino masses are always universal in this superstring 
construction (see Eq. (j2HJ)), this is not in fact completely true. It is true that they 
are universal at tree level because in this case the gauge kinetic function in Eq. (|30|) 
is always proportional to the dilaton. However, threshold corrections turn out to be 
crucial in the sin 9 — > limit, when SUSY breaking is not dominated by the dilaton. 
In this case one can observe that gaugino masses turn out to be non-universal, and 
a different phenomenology is obtained. The relic neutralino density was analyzed 
in this scenario in Ref. 230. More recently, the relic density was also analyzed in 
the general context when the soft Lagrangian is dominated by loop contributions 
in Ref. 12311 For another recent work studying also direct and indirect detection in 
these scenarios see Ref. 12321 

Let us finally remark that an analysis of the dark matter using a specific mecha- 
nism for the breaking of SUSY has also been carried out in the literature. In partic- 
ular, it was pointed out in Ref. 233 that, since scalar masses are much larger than 
gaugino masses when SUSY is spontaneously broken by non-perturbative hidden- 
sector gaugino condensates, the relic abundance of the neutralinos is too large and 
incompatible with the astrophysical observations in most of the parameter space. 

Summarizing, in the specific compactification models of the hctcrotic superstring 
studied here, although the soft terms can be non-universal, the final cross section 
is generically small and similar to the one of the mSUGRA scenario. 



5.3.2. D-brane constructions 

D-brane constructions are explicit scenarios where two of the interesting situations 
studied in Section 5.2, non-universality and intermediate scales, may occur. Con- 
cerning the latter, it was recently rea lized that the string scale may be anywhere 
between the weak and the Plank scale 1223 For instance, embedding the SM inside 
D3-branes in type I strings, the string scale is given by 

Mj = aMp !! ncfc M e 3 , (38) 
V 2 

where a is the gauge coupling and M c is the compactification scale. Thus one gets 
Mi « 10 10 - 12 GeV with M c w lO 8 " 10 GeV. 

As mentioned in Section 5. 2.2, t hese intermediate scales are interesting from 
the phenomenological viewpoint I I They are also interesting from the theoretical 
viewpoint. For example, in supergravity models supersymmetry can be sponta- 
neously broken in a hidden sector of the theory and the gravitino mass, which sets 
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the overall scale of the soft terms, is given by: 

F 

m 3/2 » jz , (39) 

JVlpianck 

where F is the auxiliary field whose vacuum expectation value breaks super- 
symmetry. Since in supergravity one would expect F w Mp lanck , one obtains 
m 3/2 ~ M Planck and therefore the hierarchy problem solved in principle by super- 
symmetry would be re-introduced, unless non-perturbative effects such as gaugino 
condensation produce F « M\nrMpi anc i-. However, if the scale of the fundamen- 
tal theory is Mj sa 10 11 GeV instead of Mpi anc k, then F ss Mj and one gets 
m 3/2 ~ Afw m a natural way, without invoking any hierarchically suppressed 
non-perturbative effect. In the above example with a D3-brane, with a modest 
input hierarchy w between string and compactification scales, Mi w 10 11 GeV and 
M c w 10 9 GeV, one obtains the desired hierarchy Myy /Mpi anc k ~ 10~ 16 . 

The first attempts to study dark matter within these constructions were car- 
Ile d out in sc enarios with the unification scale Mqut ~ 10 16 GeV as the initial 
scale EZE0E2 anc l dilaton-dominated SUSY- breaking scenarios with an intermedi- 
ate scale as the initial scale ESI. However, the important issue of the D-brane origin 
of the U(l)y gauge group as a combination of other J7(l)'s and its influence on the 
matter distribution in these scenarios was not included in those analyses. When this 
is taken into account, interesting results are obtained M * particular, scenarios 
with the gauge group and particle content of the SUSY standard model lead natu- 
rally to intermediate values for the string scale, in order to reproduce the value of 
gauge couplings deduced from experiments. In addition, the soft terms turn out to 
be generically non universal. Due to these results, in principle, large cross sections 
can be obtained. 

Let us first recall that there are two possible avenues to construct the supersym- 
metric standard model with D-branes: (a) The SU(3) C , SU(2)l and U(l)y groups 
of the standard model come from different sets of Dp-branes. (b) They come from 
the same set of Dp-branes. In the first scenario, it is worth remarking the difficulty 
of obtaining three copies of quarks and leptons if the gauge groups are attached to 
different sets of Dp-branes Thus whether or not this scenario may arise from dif- 
ferent sets of Dp-branes in explicit string constructions is an important issue which 
is worth attacking in the future. Concerning the other scenario (b), models with 
the gauge group of the standard model and three families of particles have been 
explicitly built. Let us concentrate first on scena rio fa). 

Consider for example a type I string scenario ^2 wnere the gauge group U(3) x 
U(2) x i7(l), giving rise to SU (3) x SU(2) x U(l) 3 , arises from three different types of 
D-branes as shown schematically in Fig. 1251 where open strings starting and ending 

w It is worth noticing, however, that those values would imply Re(S') = 1/ct ~ 24 and Re(T) = 
a ("SF") ~ W®, i-e. one has again a hierarchy problem but now for the VEVs of the fields that 
one has to determine dynamically. 
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Q „ 




U(3) c U(1) U(2) L 

Fig. 25. A generic D-brane scenario giving rise to the gauge bosons and matter of the standard 
model. It contains three Dp3-branes, two Dp2-branes and one Dpi-brane, where pjv may be either 
9 and 5i or 3 and 7j. The presence of extra D-branes, say Dq-branes, is also necessary. For each 
set the Dpjv-branes are in fact on the top of each other. 



on the same sets of Dp^r-branes give rise to the gauge bosons of the standard model. 
For the sake of visualization each set is depicted at parallel locations, but in fact 
they are intersecting each other. Let us recall that the presence of extra D-branes, 
say Dg-branes, is also necessary as explained in Ref. ^2 

Since the standard model gauge group arises from different types of D-branes, 
the gauge kinetic functions associated with them are different in general, and there- 
fore the gauge couplings are non-universal. On the other hand, f/(l)y is a linear 
combination of the three U(l) gauge groups arising from £7(3), U(2) and U(l) 
within the three different D-branes (with the extra Z7(l)'s anomalous and with the 
associated gauge bosons with masses of the order of the string scale Mr). This 
implies 



a Y {Mr) ai(M r ) a a (Mj) 3a 3 (M / ) ' v ; 

where o>k correspond to the gauge couplings of the U(k) branes. As shown in R ef. 1411 
in order to reproduce the low-energy value of the gauge couplings deduced from 
experiments, the above equation leads to solutions with the string scale Mj ss 
10 io-i2 GeV _ Thif , scenario is shown in Fig. ESI for Mi = 10 12 GeV. Notice that 
it is similar to the one shown in Fig. I17h . when discussing intermediate scales in 
mSUGRA, but with the qualitatively difference that the running of the U(l)y gauge 
coupling must fulfil relation H40|) . 

The analysis of the soft terms has been done under the assumption that only 
the dilaton (S) and moduli (Tj) fields contribute to SUSY breaking and it has 
been found that these soft terms are generically non-universal. Considering the 
assignment of gauge bosons and matter of Fig. 1251 and using the standard parame- 
terization of Eq. (|32l) . one is able to obtain the following soft terms^5 The gaugino 
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Fig. 26. Running of the gauge couplings of the MSSM with energy Q embedding the gauge 
groups within different sets of Dp-brancs (solid lines). For comparison the running of the MSSM 
couplings with the usual normalization factor for the hypercharge, 3/5, is also shown with dashed 
lines. An overall SUSY scale of 1 TeV is used. 



masses associated with the three gauge groups of the standard model are given by 

M 3 = V3m 3 / 2 sin0 , 
M 2 = VSm 3 /2 ©i cos0 , 



, , - , . /2 © 3 cos ©icos0 2 sin0 
M Y = V3m 3/2 a y (M/) [ + — ^— T + 



ai(Mj) a 2 (M/) 3a 3 (M/) 
The soft scalar masses of the three families are given by 

3 



m Q L 



m d R 



m„ „ = rn 



ml/2 



3/2 
2 



3/2 



m 



m 



3/2 



l 3/2 



m 



" l 3/2 



1 - Qf) cos 2 

el) cos 2 e 
e\) cos 2 e 

sin 2 6 + ej cos 2 (9) 
sin 2 9 + ej cos 2 0) 
sin 2 6 + 2 cos 2 0) 



(41) 



(42) 
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Fig. 27. 



Scatter plot of the scalar neutralino-proton cross section <r-o_ 



neutralino mass m-o in the D-brane scenario with the string scale Mj 

*1 



as a function of the 
10 12 GeV discussed 

in the text, and for tan/3 = 10 and 15. Only the big (red and blue) dots fulfil b — > and 
</ M — 2 constraints. The red ones correspond to points with > 114 GeV whereas the blue ones 
correspond to points with 91 < < 114 GeV. DAMA and CDMS current experimental limits 
and projected GENIUS limits are also shown. 



where e.g. ur denotes the three family squarks ur 7 cr, tji. Finally the trilinear 
parameters of the three families are 

v/3 

A u = —m 3/2 [(0 2 - 0i - 6 3 ) cos 9 ~ sine? ] , 



A d 



v/3 







™3/2 [(@3 - ©1 ~ ©2) cosf? - sinf? ] 



(43) 



Although these formulas for the soft terms imply that one has in principle five 
free parameters, m 3 / 2 , and &i with i = 1,2, 3, due to relation | | 2 = 1 only 
four of them are independent. In the analysis the parameters 9 and 0i are varied 
in the whole allowed range, < 9 < 2tt, — 1 < Oj < 1. For the gravitino mass, 
m 3/2 < 300 GeV is taken. Concerning Yukawa couplings, their values are fixed 
imposing the correct fermion mass spectrum at low energies, i.e., one is assuming 
that Yukawa structures of D-brane scenarios give rise to those values. 

Fig. 1771 displays a scatter plot of function of the neutralino mass m^o 

for a scanning of the parameter space discussed above ^0. Two different values of 
tan j3, 10 and 15, are shown. Although the astrophysical and UFB constraints have 
not been taken into account in the analysis, and sin 2 (a — 0) has not been computed 
explicitly, several conclusions can already be drawn. Regions of the parameter space 
consistent with DAMA limits exist, but the b — > sj and — 2 constraints forbid 
most of them. The latter are shown with small (green) points, and they have 91 < 
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Fig. 28. Schematic of heterotic M-Theory. 



in h < 114 GcV. In Fig. [57| only regions with big (red and blue) dots fulfil the above 
mentioned constraints. The red ones correspond to points with m/, > 114 GeV 
whereas the blue ones correspond to points with 91 < nth < 114 GeV. It is worth 
noticing that the larger tan/3 is, the smaller the regions allowed by the experimental 
constraints become. For example, increasing tan/? the value of — 2 turns out to 
be larger and may exceed the experimental bounds x . 

Let us finally mention that other examples with the standard model gauge group 
embedded in D-b rane s in a different way, or with larger values of the string scale, 
can also be found The dark matter cross section is qualitatively similar. 

Concerning scenario (b) 1 where all gauge groups of the standard model are 
embedded within the same set of D-branes, and therefore with gauge coupling 
unification, it can be analyzed similarly to scenario (a). However, in this case the 
value of the cross section is generically smaller. As mentioned above, in the context 
of mSUGRA with an intermediate scale, this is due to the different values of a's at 
the string scale in both types of scenarios. 

Summarizing, in D-brane configurations from type I string non-universal soft 
terms and intermediate scales arise naturally. Although we saw in Subsection 5.2 
that for some values of the parameters this may be interesting in order to increase 
the value of the cross section, in our case these values are in such a way that the 
cross section turns out to be essentially below DAMA. However, it is worth noticing 
that regions accessible for future experiments are present. 

5.4. M-Theory predictions for the neutralino-nucleon cross section 

The proposal of 11-dimensional M-Theory as a fundamental theory which contains 
the five 10-dimensional superstring theories, as well as 11-dimensional SUGRA, as 



x However, for larger values of the gravitino mass, in some special regions of the parameter space 
the allowed points may even enter in DAMA 1222) 
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different vac ua of its moduli space has motivated a large amount of phenomenolog- 
ical analyses ESI Xhe cornerstone of most of these works is the construction due to 
Hofava and Witten, who showed that the low-energy limit of M- Theory, compact- 
ified on a S 1 /Z2 orbifold with Eg gauge multiplets on each of the 10-dimensional 
orbifold fixed planes, was indeed the strong coupling limit of the Eg x Eg heterotic 
string theory. This is shown schematically in Fig. 1281 

A resulting 4-dimensional SUGRA can be obtained if six of the remaining di- 
mensions are compactified in a Calabi-Yau manifold. A certain number of virtues 
were noticed in this theory. The most relevant one was the possibility of tuning 
the 1 1-dimensional Planck scale and the orbifold radius so that the GUT-scale, 
Mgut ~ 2 x 10 16 GeV, which is here identified with the inverse of the Calabi-Yau 
volume, and the Planck scale, Mpi anc k — 1.2 x 10 19 GeV, were recovered. 

On the other hand, th e str ucture of the soft SUSY-breaking terms has been 
determined with the result E2S1 

V3m 3/2 / 1 
M = — sin 6* H — p£o cos 



1 + eo V VS 

to (6 + e Q ) sin 2 9 + (3 + 2e ) cos 2 9 - 2\/3eo sin 9 cos 9 



2 2 6m 3/2 

m = m 



3/2 (3 + £o ) 2 



V3: 



111 



3/2 



(3-2e o )sin6» + \/3eoC0s6'j , (44) 



where — 1 < eo < 1, and only one modulus T (also valid in the overall modulus 
case), with the standard parameterization F s = y/3(S + S*)m 3 / 2 sm F T — 
(T+T*)rn 3 / 2 cos 9, has been assumed. Let us remark that this assumption about the 
Calabi-Yau compactification leads to interesting phenomenological virtues EMI l n 
particular, the soft SUSY-breaking terms are automatically universal, and therefore 
the presence of dangerous FCNC is avoided. Examples of such compactifications 
exist, as e.g. the quintic hypersurface CP 4 . Although these spaces were also known 
in the context of the weakly-coupled heterotic string, the novel fact in heterotic 
M-thcory is that model building is relatively simple, and the construction of three 
generation models might be considerably easy. 

Using the above results for the soft terms, one can analyze how compatible is 
the parameter space of heterotic M-theory with the sensitivity of current dark 
matter detectors. Several analyses of dark matter in M-theory were performed 
| 23^ | 23^ | 2() | 2^7 I p a y m g special attention to the calculation of the relic density. Since 
in this scenario the soft terms (|44|l are universal, and it is very unnatural to obtain 
low scales 1^21 all examples concerning dark matter can be considered as a subset 
into the parameter space of mSUGRA with a GUT scale. In this sense, one should 
not expect to obtain high values for er^o_ p . Let us review this result. 

The soft terms are expressed in terms of three free parameters: the gravitino 
mass TO3/2, the Goldstino angle, 9, and the parameter eo- The values of the 
eo parameter will be chosen in order to guarantee the correct GUT scale (i.e., 
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Fig. 29. a) Neutralino-nucleon cross-section versus neutralino mass for —0.6 < eo < —0.1. b) 
The same for 0.1 < eo < 1. In both cases only the big (black and dark gray) dots fulfill both 
b — ► S7 and cj M — 2 constraints. The black ones correspond to points with > 114 GeV, whereas 
the dark gray ones correspond to points with 91 < < 114 GeV. Current DAMA and CDMS 
limits and the projected GENIUS limit are shown. 

—0.6 < eo < —0.1 and 0.1 < eo < 1 for non-standard and standard embeddings, 
respectively). We will take 7713/2 = 300 GeV, and tan/3 = 10, performing a vari- 
ation of the goldstino angle, 8, in [0, 27r). Both cases are depicted in Fig. 1291 
Although the astrophysical and UFB constraints have not been taken into account 
in this analysis, several conclusions can already be drawn. The experimental con- 
straints, b — > S7 and g^ — 2, put severe bounds, but still neutralinos as light as 
~ 100 GeV can be obtained. Once the lower limit on the Higgs mass is applied, the 
cross-section is as small as cr%o_ p ~ 3 x 10 -9 pb, far beyond the reach of present 
detectors, and close to the lower limit of the projected GENIUS. Let us recall that 
all the points represented satisfy the experimental constraints on the lower masses 
of the supersymmetric particles and satisfy rrih > 91 GeV. Small (light gray) dots 
represent points not fulfilling the b — > sj constraint. Large dots do satisfy that con- 
straint, and among these, dark gray points have 91 GeV< mh < 114 GeV, while 
black dots satisfy the stronger lower bound for the Higgs mass > 114 GeV. Of 
course, as discussed above, this constraint on the Higgs mass holds in general for 
the cases with universal soft terms for tan (3 < 50 and therefore it is the one we 
should consider here. However, due to the strong restrictions it imposes it is shown 
explicitly. 

Although the predicted values for the cross-section increase, in principle, when 
larger values of tan (3 are taken into account, the experimental bounds become much 
more important in these cases (especially those corresponding to b — > S7 and g^ — 2), 
excluding larger regions in the parameter space and thus forbidding large values of 
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Finally, it is worth noticing that non-perturbative objects of M-Theory, such as 
M5-branes, can be shown to survive the orbifold projection of Hofava-Witten con- 
struction under certain circumstances, permitting much more freedom to play with 
gauge groups and with the matter fields that appear. In this context, dark matter 
implications of vacua with five-branes were investigated in the limit where the five- 
branc modulus, Z, is the only one responsible for the breaking of supersymmetry 
^21. However, the authors used prev ious soft-terms computed in the literature, where 
some corrections were not included* ' ■ " ■ 'l The soft terms are now more involved than 
those in Eq. I|44|) and include the F terms associated with S, T and Z. A similar 
analysis as above, concerning dark matter detection, can be carried out in this case 

The highest value of the cross section, cr%o_ p ~ 10~ 8 pb, can be obtained in the 
special case in which the five-brane modulus is the only one responsible for SUSY 
breaking. 

Summarizing, as it could be expected from the universality of the soft terms, 
the predictions of heterotic M-theory with one modulus for the neutralino-nucleon 
cross section are too low to be probed by the present dark matter detectors. Only 
future experiments, as e.g. GENIUS, would be able to explore such low values. 



6. Conclusions 

Nowadays there is overwhelming evidence that most of the mass in the universe 
(90% and probably more) is some (unknown) non-luminous 'dark matter'. At galac- 
tic and cosmological scales it only manifests through its gravitational interactions 
with ordinary matter. However, at microscopical scales it might manifest through 
weak interactions, and this raises the hope that it may be detected in low-energy 
particle physics experiments. 

Plausible candidates for dark matter are Weakly Interacting Massive Parti- 
cles, the so-called WIMPs. They are very interesting because they can be present 
in the right amount to explain the observed matter density of the Universe 
0.1 < ilh 2 < 0.3. The leading candidate for WIMP is the so-called neutralino, a 
particle predicted by the supersymmetric extension of the standard model. These 
neutralinos are stable and therefore may be left over from the Big Bang. Thus they 
will cluster gravitationally with ordinary stars in the galactic halos, and in particu- 
lar they will be present in our own galaxy, the Milky Way. As a consequence there 
will be a flux of these dark matter particles on the Earth. 

Many underground experiments have been carried out around the world in order 
to detect this flux. One of them, the DAMA collaboration, even claims to have de- 
tected it. They obtain that the preferred range of the WIMP-nucleon cross section 
is pb 10~ 6 - 10~ 5 pb for a WIMP mass between 30 and 270 GeV. Unfortunately, 
this result is controversial because of the negative search result obtained by other 
experiments like CDMS, EDELWEISS AND ZEPLIN in the same range of param- 
eters. Thus we will have to wait for the next generation of experiments, which are 
already starting operations or in project, to obtain more information about whether 
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or not ncutralinos, or generically WIMPs, are the evasive dark matter filling the 
whole Universe. The most sensitive detector, GENIUS, will be able to test a WIMP- 
nucleon cross section as low as w 10~ 9 pb. Indeed such a sensitivity covers a large 
range of the parameter space of SUSY models with neutralinos as dark matter. 

Concerning this point, we have reviewed the known SUSY scenarios, and in 
particular how big the cross section for the direct detection of neutralinos can be. 
This analysis is crucial in order to know the possibility of detecting dark matter in 
the experiments. In particular, we have concentrated on SUGRA and superstring 
and M-Theory scenarios. 

Let us recall that the analysis has been carried out imposing the most recent 
experimental and astrophysical constraints on the parameter space. Concerning the 
former, the lower bound on the Higgs mass, the b — ► S7 branching ratio, and the 
muon g — 2 have been considered. The astrophysical bounds on the matter density 
mentioned above have also been imposed on the theoretical computation of the 
relic neutralino density, assuming thermal production. In addition, the constraints 
that the absence of dangerous charge and colour breaking minima imposes on the 
parameter space has also been taken into account. 

In the usual mSUGRA scenario, where the soft terms are assumed to be uni- 
versal, and the GUT scale is considered, the cross section is constrained to be 
&Xi-p ~ ^ x 10 -8 pb- Obviously, in this case, present experiments are not suffi- 
cient and more sensitive detectors producing further data are needed. A similar 
conclusion is obtained when an intermediate scale is considered. Although the cross 
section increases significantly, the experimental bounds impose er^o_ p < 4 x 10~ 7 
pb. And, in fact, at the end of the day, the preferred astrophysical range for the 
relic neutralino density, 0.1 < tt^oh 2 < 0.3, imposes cr^o_ p < 10~ 7 pb. Still present 
experiments are not sufficient. 

When non-universal scalars are allowed in SUGRA, for some special choices of 
the non-universality, the cross section can be increased a lot with respect to the 
universal scenario. It is even possible, for some particular values of the parame- 
ters, to find points allowed by all experimental and astrophysical constraints with 
<T^o_ p rj 10~ 6 pb, and therefore inside the DAMA area. This is similar to what 
occurs in the so-called cffMSSM scenario. For non-universal gauginos, although the 
cross section increases, the experimental bounds exclude this possibility. 

On the other hand, the low-energy limit of superstring theory and M-theory 
is SUGRA, and therefore the neutralino will also be a candidate for dark mat- 
ter in these scenarios. In the context of superstring theory we have reviewed two 
interesting constructions, the perturbative heterotic superstring and D-branes con- 
figurations from type I string, where the soft terms can be computed explicitly in 
some models. In the former, although the soft terms can be non- universal, the final 
cross section is similar to the one of the mSUGRA scenario. In the latter, in ad- 
dition to non-universality, intermediate scales arise naturally, but again the cross 
section is in general small (although some regions may even be compatible with 
DAMA). 
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Finally, due to the universality of the soft SUSY-breaking terms in the heterotic 
M-Theory scenario analyzed with only one modulus, and the fact that the most 
natural value for the initial scale is of order 10 16 GeV, the parameter space can be 
considered as a subset of mSUGRA. Therefore, the predicted cross-section is very 
low, cr^o_ p < 10 -8 pb, far beyond the reach of the present dark matter experiments. 

In summary, underground physics as the one discussed here is crucial in order 
to detect dark matter. Even if ncutralinos are discovered first at future particle 
accelerators such as LHC, only their direct detection due to their presence in our 
galactic halo will confirm that they are the sought-after dark matter of the Uni- 
verse. For that, many new underground experiments are already starting operations 
or in project, and they will be able to cover an important range of the parameter 
space of SUSY models. 
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